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Abstract: 

We consider two-loop leading and next-to-leading logarithmic virtual corrections to 
arbitrary processes with external massless fermions in the electroweak Standard Model 
at energies well above the electroweak scale. Using the sector-decomposition method 
and alternatively the strategy of regions we calculate the mass singularities that arise 
as logarithms of Q 2 /M^, where Q is the energy scale of the considered process, and \/e 
poles in D = 4 — 2e dimensions, to one- and two-loop next-to- leading logarithmic accuracy. 
The derivations are performed within the complete electroweak theory with spontaneous 
symmetry breaking. Our results indicate a close analogy between the form of two-loop 
electroweak logarithmic corrections and the singular structure of scattering amplitudes 
in massless QCD. We find agreement with the resummation prescriptions that have been 
proposed in the literature based on a symmetric SU(2) x U(l) theory matched with QED 
at the electroweak scale and provide new next-to-leading contributions proportional to 
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1 Introduction 



The electroweak radiative corrections to high-energy processes are characterized by 
the presence of logarithms of the type ln(Q 2 /M 2 ), which involve the ratio of the typical 
scattering energy Q over the gauge-boson mass scale M = M w ~ M z [1] . These logarith- 
mic corrections affect every reaction that involves electroweakly interacting particles and 
has a characteristic scale Q 3> M. They start to be sizable at energies of the order of a few 
hundred GeV and their impact increases with energy. At the LHC [2] and future lepton 
colliders [3-5], where scattering energies of the order of 1 TeV will be reached, logarithmic 
electroweak effects can amount to tens of per cent at one loop and several per cent at 
two loops. Thus, this class of corrections will be very important for the interpretation of 
precision measurements at future high-energy colliders. 

For sufficiently high Q, terms of 0(M 2 /Q 2 ) become negligible and, at I loops, the 
electroweak corrections assume the form of a tower of logarithms, 



where the leading logarithms (LLs), also known as Sudakov logarithms [6], have power 
j = 21, and the subleading terms with j = 21 — 1, 21 — 2, . . . are denoted as next-to- 
leading logarithms (NLLs), next-to- next-to-leading logarithms (NNLLs), and so on. The 
complete asymptotic limit includes all logarithmic contributions (j > 0) as well as terms 
that are not logarithmically enhanced at high energy (j = 0)0 

Electroweak logarithmic corrections have a twofold origin. On the one hand they can 
appear as terms of the form ln(Q 2 / /i^) resulting from the renormalization of ultraviolet 
(UV) singularities at the renormalization scale /xr ~ M. These logarithms can easily be 
controlled through the running of the coupling constants. The other source of electroweak 
logarithms are mass singularities, i.e. logarithms of the form \n(Q 2 /M 2 ) that are formally 
singular in the limit where the gauge-boson masses tend to zero. As is well known, in 
gauge theories, mass singularities result from the interactions of the initial- and final- 
state particles with soft and/or collinear gauge bosons. This permits, in principle, to 
treat mass singularities in a process-independent way and to derive universal properties 
of mass-singular logarithmic corrections. 

At one loop, it was proven that the electroweak LLs and NLLs are universal, and a gen- 
eral formula was derived that expresses the logarithmic corrections to arbitrary processes 
in terms of the electroweak quantum numbers of the initial- and final-state particles [7,8]. 
In recent years, the impact of one-loop electroweak corrections was studied in detail for 
various specific processes at high-energy colliders [9-15]. At the LHC, in general, large 
negative corrections are observed that appear at transverse momenta around 100 GeV 
and grow with px- Depending on the process, the size of these corrections can reach 
up to 10-40% at p T ~ 500-1000 GeV. In Refs. [12,14], the predictions based on NLL 
and NNLL high-energy approximations were compared with exact one-loop calculations. 

1 In the literature, such terms are sometimes denoted as "constants" since they do not involve log- 
arithms that grow with energy. However, in general they are functions of the ratios of kinematical 
invariants, which depend on the scattering angles. At one loop such terms are formally classified as 




with 




(1.1) 



NNLLs. 
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In both cases it turned out that in the high-pp region, where the corrections are large, 
the NNLL predictions deviate from the exact calculation by much less than 1%, i.e. the 
numerical effect of 0(M 2 / Q 2 ) contributions is negligible. Also the NLL approximation 
provides a correct description of the bulk of the corrections and their energy dependence. 
However, it was found that the actual precision of the NLL approximation, better than 
1% for pp — » jet + Z/7 [14] and about 5% for pp — > W7 [12], depends relatively strongly 
on the process. As it was shown for pp — ► jet + Z/7, also the two-loop electroweak loga- 
rithms can have an impact at the several per-cent level on high-px measurements at the 
LHC [14]. 

In the recent years, the properties of electroweak logarithmic corrections beyond one 
loop have been studied with two complementary approaches^] On the one hand, evolution 
equations, which are well known in QED and QCD, have been applied to the electroweak 
theory in order to obtain the higher-order terms through a resummation of the one-loop 
logarithm^ [17-23]. O n the other hand, explicit diagrammatic calculations based on the 
electroweak Feynman rules have been performed [24-31]. 

Fadin et al. [17] have resummed the electroweak LL corrections to arbitrary matrix ele- 
ments by means of the infrared evolution equation (IREE), and this approach has been ex- 
tended by Melles to the NLL terms for arbitrary processes [18-21]. In Refs. [22,23] Kiihn 
et al. have resummed the logarithmic corrections to neutral-current massless 4-fermion 
processes up to the NNLL level. These resummations [17-23] are based on the IREE [17], 
which describes the all-order LL dependence of matrix elements with respect to the cut- 
off parameter fi±. This cut-off fixes the minimal transverse momentum for the gauge 
bosons that couple to the initial- and final-state particles and acts as a regulator of mass 
singularities. The IREE, which was originally derived within symmetric gauge theories 
(QED and QCD), has been applied to the spontaneously broken electroweak theory under 
the assumption that this latter can be split into two regimes with exact gauge symmetry. 
In practice, in the regime Q > fi± > M the masses of the electroweak gauge bosons, 
which result from the breaking of gauge symmetry, are supposed to be negligible and 
exact SU(2)xU(l) symmetry is assumed. Instead, in the regime M > fi± the weak gauge 
bosons are supposed to be frozen owing to their masses such that only photons contribute 
to the evolution, and this is characterized by exact U(l) em symmetry. 

As a result, the resummed electroweak corrections factorize into two parts correspond- 
ing to these two regimes: (i) a symmetric-electroweak part that can be computed within 
a symmetric SU(2)xU(l) gauge theory using the same cut-off parameter M to regularize 
the mass singularities that result from all gauge bosons, i.e. assuming that the photon 
is as heavy as the weak gauge bosons, and (ii) an electromagnetic part that originates 
from the mass gap between photons and massive gauge bosons and can be computed 
within QED. The electromagnetic part contains divergences that are due to massless 
photons and depend on the scheme adopted for their regularization. These divergences 
cancel when the contributions of virtual and real photons are combined, and then the 

2 For recent developments in the exact numerical calculation of complete two- loop integrals see, for 
instance, Ref. [16]. 

3 In order to resum the LLs and NLLs it is sufficient to determine the kernel of the evolution equations 
to one-loop accuracy. However, starting from the NNLLs also the two-loop contributions to the /3-function 
and to the anomalous dimensions are needed. 
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electromagnetic contribution depends on the cuts imposed on real photons. Instead the 
symmetric-electroweak part involves logarithms of the form ln(Q 2 /M 2 ) which are for- 
mally mass singular but numerically finite since M does not vanish. These logarithms 
are present in all physical observables that are exclusive with respect to real radiation of 
Z and W bosons. Moreover, the electroweak logarithms remain present even in inclusive 
observables due to a lack of cancellation between virtual and real contributions from elec- 
troweak gauge bosons. This is due to the fact that the conditions of the Bloch-Nordsieck 
theorem are not fulfilled since fermions and gauge bosons carry non-abelian weak-isospin 
charges [32]. 

The resummations [17-23] rely on the assumption that all relevant implications of 
electroweak symmetry breaking are correctly taken into account by simply splitting the 
evolution equation into two regimes with exact gauge symmetry. In particular, the follow- 
ing assumptions are explicitly or implicitly made, (i) In the massless limit M/Q — > 0, all 
couplings with mass dimension, which originate from symmetry breaking, are neglected, 
(ii) The weak-boson masses are introduced in the corresponding propagators as regulators 
of soft and collinear singularities from W and Z bosons without spontaneous symmetry 
breaking. Since these masses are of the same order, one uses M w = M z as an approxima- 
tion, (iii) The regimes above and below the electroweak scale are treated as an unmixed 
SU(2)xU(l) theory and QED, respectively, and mixing effects in the gauge sector, i.e. 
the interaction of photons with W bosons, are neglected. 

It is important to understand to which extent the above assumptions are legitimate 
and whether the resulting resummation prescriptions are correct. This can be achieved by 
explicit diagrammatic two-loop calculations based on the electroweak Lagrangian, where 
all effects related to spontaneous symmetry breaking are consistently taken into account. 
At the LL level, these checks have already been completed. A calculation of the massless 
fermionic singlet form factor [24, 25] and then a Coulomb-gauge calculation for arbi- 
trary processes [26] have supported the exponentiation of the electroweak LLs predicted 
by the IREE. Also the angular-dependent subset of the NLL corrections for arbitrary 
processes [27] has been shown to be consistent with the exponentiation anticipated in 
Refs. [21-23]. At present the complete set of electroweak NLLs is known only for the 
so-called fermionic form factor, which corresponds to the gluon-fermion-antifermion ver- 
tex [28]. This calculation has confirmed that the IREE approach provides a correct 
description of the terms that do not depend on the difference between the masses of 
the Z and the W bosons and has provided first insights into the behaviour of the NLL 
terms of the type a 2 ln(M|/M w ) ln 3 (Q 2 /M 2 ). The complete tower of two-loop logarithms 
for the fermionic form factor has been computed in Refs. [29,30] adopting an unmixed 
SU(2)xU(l) theory with Mz = Mw as approximation. Within this framework it was 
found that the soft-collinear singularities resulting from massless photons factorize as 
suggested by the IREE and that, in absence of mixing, symmetry breaking is negligible 
up to the NNLL level and its first non-trivial effects appear through a Higgs-mass depen- 
dence of the NNNLL terms. Using QCD resummation techniques, the NNNLL two-loop 
results for the fermionic form-factor have been extended to the amplitude for neutral- 
current four-fermion scattering [30,31]. Here terms of the type In (M|/M w ) were included 
through an expansion in = 1 — M w /M| up to the first order in s^. 



3 



In this paper we develop a formalism to derive virtual electroweak two-loop LL and 
NLL corrections for arbitrary processes and apply it to the case of massless n-fermion 
reactions. The calculation is performed diagrammatically using the electroweak Feynman 
rules in the 't Hooft-Feynman gauge. The mass-singular logarithms are extracted from the 
relevant Feynman diagrams by means of a soft-collinear approximation for the interaction 
of initial- and final-state fermions with gauge bosons. Using Ward identities we prove that 
the LL and NLL mass singularities for generic n-fermion amplitudes factorize from the 
corresponding Born amplitudes. All relevant loop integrals are evaluated in the high- 
energy region Q ^> M to NLL accuracy using an automatized algorithm based on the 
sector-decomposition technique [33] and alternatively the method of expansion by regions 
combined with Mellin-Barnes representations [34]. We do not assume Mi = Mw and we 
include all relevant contributions depending on the difference Mz — Mw- In addition to 
the logarithms of the type (11.11) . which arise from massive virtual particles, we include also 
mass singularities from massless virtual photons. The latter are regularized dimensionally 
and arise as 1/e poles in D = 4 — 2e dimensions. In our result the photonic singularities 
are factorized in a gauge-invariant electromagnetic term. The remaining part of the 
corrections, which is also gauge invariant and does not depend on the scheme adopted to 
regularize photonic singularities, contains only finite log(Q 2 /M 2 ) terms. The divergences 
contained in the electromagnetic term cancel if real-photon emission is included. 

The paper is organized as follows. Section [2] contains definitions and conventions 
used in the calculation. In Sect. [3] we identify the diagrams that produce ultraviolet and 
mass singularities, split them into factorizable and non-factorizable parts and discuss our 
method to treat these contributions. Using collinear Ward identities, we prove in Sect. H] 
that all non-factorizable contributions cancel. Explicit expressions for the factorizable 
contributions and for the counterterms needed for the renormalization are provided in 
Sects. [5] and O respectively. The complete one- and two-loop results are presented in 
Sect. [7] and discussed in Sect. |HJ Our conclusions are contained in Sect. and various 
appendices are devoted to technical aspects of the calculation. 



2 Definitions and conventions 



We consider a generic n — > process involving an even number n of polarized fermionic 
particles, 

(flipi) . ..<fn(Pn) -> 0. (2.1) 

The symbols (fi represent n/2 antifermions and n/2 fermions: (fi = f£j for i = 1, . . . , n/2 
and ifi = f£* for i — n/2 + 1, . . . , n. These particles are massless chiral eigenstates with 



m k = and chirality «j = R or L. In practice every fermion or antifermion can be a 



lepton or a light quark. The indices <jj characterize the lepton/quark nature, the isospin 
and the generation of the fermions, j ai = u e , e, u,d, . . .. 
The matrix element for the process (12.11) reads 



M 



n/2 



i=l 



G£l-£n(p h ...,p n ) 



n u (pv k j) 



j=n/2+l 



(2.2) 
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where G-i' "-« is the corresponding truncated Green function. The spinors for chiral 
fermions and antifermions fulfil 

u p u{p, re) = S Kp u(p, re), v{p, k)u p = 5 Kp v(p, re) (2.3) 

for p, re = R, L, and 

= = ^(l+7 5 ), u L = u R = ^(l-7 5 )- (2.4) 

The matrix elements ( 12.2ft are often abbreviated as .M = A^^ 1 '"^™. 

The amplitudes for physical scattering processes, i.e. 2 — > n — 2 reactions, are easily 
obtained from our results for n — > reactions using crossing symmetry. 

2.1 Perturbative and asymptotic expansions 

For the perturbative expansion of the matrix elements in a = e 2 /(4n), where e is the 
electromagnetic coupling constant, we write 

oo 

M = J2 M n (2-5) 

1=0 

with 



and 



M, = ( ^ ) A1„ (2.6) 



Note that, for convenience, in D = 4 — 2e dimensions we include in the definition of a e 
a normalization factor depending on e, the scale /xd of dimensional regularization, the 
characteristic energy Q of the scattering process and Euler's constant 7e- For p^ = 
Q 2 e 7E /(47r) we have a e = a. 

The electroweak corrections are evaluated in the region where all kinematical invari- 
ants, rjk = (pj + Pk) 2 , are much larger than the squared masses of the heavy particles 
that enter the loops, 

\r jk \ ~Q 2 ^M^~Ml~ml~Ml. (2.8) 

This means that we consider a situation with essentially two different scales, one fixed by 
Q, which can be identified with (the square root of the absolute value of) one of the kine- 



matical invariants, Q = y\rjk\, and one fixed by the heavy masses of the Standard Model. 
In this region, the electroweak corrections are dominated by mass-singular logarithms, 
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and logarithms of UV origin. Mass singularities that originate from soft and collinear 
massless photons and UV singularities are regularized dimensionally and give rise to 1/e 
poles in D = 4 — 2e dimensions. Therefore, the /-loop contributions to the polarized 
matrix elements can be written as an expansion in L and e, 



2/ 



-H = E E M^ n e n L m+n . (2.10) 



m=0 n- 



The logarithmic terms in f 1 2 . 1 j) are classified according to their degree of singularity, 
defined as as the total power m of logarithms L and 1/e poles. The maximal degree of 
singularity at /-loop level is m = 21 and the corresponding terms are denoted as leading 
logarithms (LLs). The terms with m = 2/ — 1, 2/ — 2, . . . represent the next-to-leading 
logarithms (NLLs), next-to-next-to- leading logarithms (NNLLs), and so on. 

In this paper we systematically neglect mass-suppressed corrections of order M^/Q 2 
and we calculate the one- and two-loop corrections to NLL accuracy, i.e. including LLs 
and NLLs. For this approximation we use the symbol = . The two-loop corrections are 
expanded in e up to the finite terms, i.e. contributions of order e°L 4 and e°L 3 . Instead, 
the one- loop corrections are expanded up to order e 2 , i.e. including terms of order e 2 L A 
and e 2 L 3 . These higher-order terms in the e-expansion must be taken into account when 
expressing two-loop mass singularities in terms of one-loop ones. 

Since the loop corrections depend on various masses, M\y ~ Mz ~ m t ~ Mh, and 
different invariants r^, the coefficients M lmn in (I2.10p depend on mass ratios and ratios 
of invariants^ Actually, as we will see in the final result, they depend only on logarithms 
of these ratios, 



''^m)' (2 - n) 



The logarithms of —r^ are extracted in the region < 0, where the corrections are 
real. The imaginary parts that arise in the physical region can be obtained via analytic 
continuation replacing by + iO. 

2.2 Symmetry breaking and mixing 

In the electroweak Standard Model, the physical gauge bosons result from the gauge 
bosons W 1 ^ 2 ^ 3 and B associated with the SU(2) and U(l) gauge groups via the 
unitary transformation 



W ± 



-±=(w^TiW 2 ), Z^ = c w Wl + s w B^ A^ = -s w Wl + c w B^ (2.12) 



where c w = cos^w, s w = sin^w, and 9 W is the weak mixing angle. The generators 
associated with the physical gauge bosons are related to the weak isospin generators T l 
and the weak hypercharge Y through 



el w± = -y= (T 1 ± iT 2 ) , el z = g 2 c w T 3 - gis w —, eI A = -g 2 s w T 3 - gic w —, 



(2.13) 



This dependence only appears at the next-to- leading level, i.e. for m = 21 — 1. 
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where g% and g 2 are the coupling constants associated with the U(l) and SU(2) groups, 
respectively. The electromagnetic charge operator is defined as 

Q = -I A , (2.14) 

and the generator associated with the Z boson can be expressed as 

el z = ^T 3 -e^Q. (2.15) 

c w c w 

The SU(2)xU(l) Casimir operator is given by 

V=A,Z,W± e v z / e 

where denotes the complex conjug ate of V, and C = Ei=i( Ti ) 2 is the Casimir oper- 
ator of the SU(2) group with eigenvalues 3/4 and for left- and right-handed fermions, 
respectively. The generators (I2.13P obey the commutation relations 



rVi T v 2 



ig 2 J2 e VlV * Va I Va , (2.17) 

V 3 =A,Z,W± 



with 



£ v lV2 v 3 = ix 



;-i) p+1 c w if v 1 v 2 v 3 = tt(zw + w~), 

(-lfs w if V 1 V 2 V 3 = ir(AW + W-), (2.18) 
otherwise, 



and (— l) p represents the sign of the permutation n. 

The SU(2) x U(l) symmetry is spontaneously broken by a scalar Higgs doublet which 
acquires a vacuum expectation value v. We parametrize the Higgs doublet in terms of 
the four degrees of freedom = H, x, 4> + , 4>~ , where <p~ = (<ft + ) + . In this representation 
the gauge-group generators are 4x4 matrices with components 1®^. ■ The gauge-boson 
masses and the mixing parameters are fixed by the condition that the gauge-boson mass 
matrix is diagonal, 

Myy, = leV {l V J V '} HH = SvvMy, (2.19) 

for V, V = A, Z, W^. The curly brackets in (12.191) denote an anticommutator. The mass 
eigenvalues are given by 

M w± = -g 2 v, M z = —g 2 v, M A = 0. (2.20) 

The weak mixing angle 8 W is related to the gauge-boson masses via 

c (221) 
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The vanishing mass of the photon is connected to the fact that the electric charge of the 
vacuum expectation value is zero. This provides the relation 

Cw</i># = s w g 2 , (2.22) 

between the weak mixing angle, the coupling constants gi, g 2 and the hypercharge Y$ of 
the Higgs doublet. Identifying e = c w g±, the Gell-Mann-Nishijima relation for general 
F$ reads Q = Y/2 + Y$T 3 . In the calculation we keep Y$ as a free parameter, which 
determines the degree of mixing in the gauge sector. In this way our analysis applies to 
the Standard Model case, Y$ = 1, as well as to the case Y$ = corresponding to an 
unmixed theory with 

s w = 0, c w = l, Z ll = W*, — 13^, M W = M Z . (2.23) 

2.3 Gauge interactions of massless fermions 

For the Feynman rules and various group-theoretical quantities we adopt the formalism 
of Ref. [8] (see Apps. A and B therein). The Feynman rules for the vector-boson-fermion- 
antifermion vertices read 




where V = A,Z,W ± , and lJ K f K are the generators that describe SU(2)xU(l) transfor- 
mations of fermions in the fundamental (n = L) or trivial (k = R) representation. For 
antifermions we have 

The chiral projectors oj k that are associated with the gauge couplings (I2.24p can easily 
be shifted along the fermionic lines using anticommutation relational until they can be 
eliminated using f)2.3p . As an example, for the coupling of a vector boson V to an incoming 
fermion ^ = f£? one obtains 

Hj '' ' i^-fJ w u(Pt, «i) with i$.<pi = ifHfiii ( 2 - 26 ) 



(Pi + q) 



where the gauge coupling depends on the chirality Ki of the fermion. Similarly, for 
the coupling of a vector boson to an incoming antifermion ipi = f£j one obtains 

- v(Pi}Ki)i e^jMi±^ ^ I\ K = I v op >= (2.27) 

5 In our calculation we use {7^, 7 s } = in D = 4 — 2e dimensions. 
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where the minus sign of the propagator is absorbed in the coupling of the antifermion. 

In our results, the matrix elements (12. 2p are often multiplied by various matrices 
resulting from the couplings of gauge bosons to the external fermions. For such expressions 
we introduce the shorthands 



mi?i? = E m^--^Q,,Q 



fk^k fkfk 1 



etc. The gauge couplings in (12.281) satisfy the commutation relations 



= ^fe E e 

V 3 =A,Z,W± 



and the e-tensor is defined in (12.181) . 

Global gauge invariance implies the charge-conservation relation 



'2.28) 



(2.29) 



■ME4 

k=l 



V 



0. 



(2.30) 



which is fulfilled up to mass-suppressed terms in the high-energy limit. 



3 Treatment of ultraviolet and mass singularities 

Large logarithms and 1/e poles originate from UV and from mass singularities. In this 
section we present the technique that we use to extract these singularities from one- and 
two-loop Feynman diagrams. In Sect. 13 .11 we identify the diagrams that are responsible 
for mass singularities within the 't Hooft-Feynman gauge and classify their contributions 
into factorizable and non-factorizable ones. In Sect. 13.21 we introduce an approximation 
that describes the exchange of virtual gauge bosons in the soft and collinear limit. Our 
treatment of UV singularities is discussed in Sect. 13.31 

3.1 Mass singularities 

In gauge theories, as is well known, mass singularities appear in loop diagrams involv- 
ing virtual gauge bosons that couple to the on-shell external legsJU These singularities 
originate from the integration over the loop momenta in the regions where the momenta 
of the virtual gauge bosons become soft and/or collinear to the external momenta. 

6 In principle it is possible to adopt particular gauge fixings in order to isolate mass singularities in 
a smaller subset of diagrams. In the Coulomb gauge, for instance, collinear singularities appear only in 
external-leg self-energy corrections [26]. However, we prefer to use the 't Hooft-Feynman gauge for our 
analysis of mass singularities. Thus, we must consider the most general set of Feynman diagrams which 
gives rise to mass singularities. 
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3.1.1 Mass-singular diagrams at one loop 

At one loop, the mass singularities of the n-fermion amplitude originate from diagrams 
of the type 

V 

* , (3.i) 

where an electroweak gauge boson, V = A,Z, W ± , couples to one of the external fermions, 
% — 1, . . . , n, and to any other line of the diagram. These diagrams can be classified into 
three types: 

(i) The diagrams where the virtual gauge boson V couples to the external leg % and to 
another external leg j with j ^ i, 




These diagrams produce single- and double-logarithmic mass singularities that orig- 
inate from the regions where the gauge-boson momentum becomes soft and/or 
collinear to the momentum of the external leg i or j. 

(ii) The external-leg self-energy insertions 

V 

• ' > • (3-3) 

These diagrams constitute a subset of the general mass-singular diagrams depicted 
in (13. ip . However, they have to be omitted in our calculation since we express 
S'-matrix elements in terms of truncated Green functions and on-shell renormalized 
fields. 

(iii) The diagrams where the virtual gauge boson V couples to the external line i and to 
an internal line, i.e. an internal propagator of the tree subdiagram that is represented 
as the grey blob in (13.1 1) . These diagrams produce only single-logarithmic mass 
singularities that originate from the region where the momenta of the gauge boson 
and the external fermion i become collinear. 

3.1.2 Factorizable and non-factorizable contributions at one loop 

The diagrams of type (i) involve an n-fermion tree subdiagram. It is thus possible 
to extract from the diagrams (I3.2p mass-singular contributions that factorize from the 
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n-fermion Born matrix element. Such contributions are called factorizable and are defined 
aa3 



• 2 A-D V- jV jV I d Q 



<Pi<Pi 



{27t) D {q 2 -M^{p i + qf(p j -qf 



X v(pj, - q)C- [l>\. ...,Pi,...,Pj,.. . , Pn )(tfi + thitUfa, Hi). 

(3.4) 

Here G- 1 "- i "'-i"'- n denotes the truncated Green function corresponding to the n-fermion 
tree subdiagram. By definition, in the factorizable contributions (F) we include only those 
parts of the above diagrams that are obtained by performing the loop integration with the 
momentum q of the gauge boson V set to zero in the tree subdiagram. This prescription 
is indicated by the label F in the tree subdiagram. In (13.41) the spinors of the external 
fermions k = 1, . . . , n with k ^ i,j are implicitly understood. 

By construction, the factorizable terms (13.41) contain all one-loop soft singularities, 
since the soft singularities originate only from the diagrams of type (13.21) in the region 
where the gauge-boson momentum tends to zero. Actually, as we will show, the fac- 
torizable contributions (13 .4p contain all one-loop mass singularities, i.e. not only all soft 
singularities but also all collinear singularities. 

The combination of all factorizable one-loop contributions is obtained by summing 
over all gauge bosons V = A, Z, W ± and external legs i, j in (13. 4p . 



= iEE E 

z i=i j=i V=A,Z,W ± 

3*A 




(3.5) 



Note that the symmetry factor 1/2 avoids double counting of the pairs of external legs 
The limit q^ — > xpf,xpj in ( 13.51) indicates that the above diagrams are evaluated 
in the approximation where the four-momentum g M of the gauge boson V is soft and/or 
collinear to one of the momenta of the external legs i and j. This approximation, which 
is defined in Sect. 13.21 permits to extract all mass singularities, which result from soft 
and/or collinear regions]^ As we will see, these terms can be expressed as products of 
n-fermion Born matrix elements and one-loop integrals. The factorizable one-loop terms 
(13. 5p are computed in Sect. 15. 11 



7 Here we consider the case where the particles i and j are a fermion and an antifermion, respectively. 
The generalization to other combinations of particles or antiparticles is obvious. 

8 Here we adopt a different approach as compared to Ref. [7]. There the diagrams of type (|3.5|) 
were evaluated using the cikonal approximation, which extracts only soft singularities, and the collinear 
singularities where treated separately. 
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The remaining one-loop mass singularities are obtained by subtracting the factorizable 
contributions (13. 4p and the external self-energy contributions (13.31) from the diagrams 

(ED, 



V 



V 



II 



M 



NF 

1 



E E 



i=l V=A,Z,W ± 



E 



n 




(3.6) 



j=i 



These contributions are called non-factorizable (NF). We note that these NF terms are 
free of soft singularities since all soft singularities are contained in the factorizable parts 
and are subtracted in ( 13. 6p . We also observe that, in contrast to ( 13.51) . the sum over 
pairs of external legs i,j appearing in front of the factorizable contributions that are 
subtracted in (13.61) is not multiplied by a symmetry factor 1/2. This is due to the fact 
that in the subtraction terms in (13. 6p we include only the contribution of the collinear 
region g M — > xpf, whereas in (13.51) every diagram contains the mass singularities resulting 
from the two regions g M — * xpf and — > xpj. 

In Sect. H]we will prove that the non-factorizable one-loop terms ( 13.61) vanish. 

3.1.3 Mass-singular diagrams at two loops 

The diagrams that give rise to NLL mass singularities at two loops, i.e. terms with 
triple and quartic logarithmic singularities, can be obtained from the one-loop diagrams 
of type (I3.2p . which produce double logarithms in the soft-collinear region, by inserting 

• a second soft and / or collinear gauge boson that couples to an external line or to the 
virtual gauge boson in (13. 2p providing an additional single or double logarithm, 

• or a self-energy subdiagram in the propagator of the virtual gauge boson in (13.21) . 
which provides an additional single logarithm. 

There are five types of such diagrams: 



Here the NLL mass singularities originate from the regions where the gauge boson V±, 
which couples to two external legs or to an external leg and another virtual gauge boson, 
is soft and collinear and the gauge bosons Vi and V3, in the first four diagrams in ( 13.71) . 
are soft and/or collinear. 




(3.7) 
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3.1.4 Factorizable and non-factorizable contributions at two loops 



The two-loop mass singularities resulting from the diagrams (13.71) are split into fac- 
torizable and non-factorizable contributions. The factorizable contributions result from 
the diagrams that contain n-fermion tree subdiagrams, 




L J / qft—txpt 1 

(3.8) 

As in the one-loop case, these factorizable contributions are defined as those parts of the 
diagrams (13.81) that are obtained by performing the loop integrations with the momenta 
of the gauge bosons Vi,V 2 ,V 3 set to zero in the tree subdiagrams. This prescription is 
indicated by the label F in the tree subdiagrams in (13.8H . 

By construction, the factorizable terms H 3 . 8 1) contain all two-loop soft-soft singulari- 
ties, i.e. the singularities that originate from the region where all gauge bosons V\, V2, V3 
are soft. Actually, as we will show, the factorizable contributions (13.81) contain all two- loop 
NLL mass singularities, i.e. not only all soft singularities but also all collinear singularities. 

The symmetry factors 1/2, 1/6 and 1/8 in (13. 8p avoid double counting in the sums 
over combinations of external legs i,j, k, I. The limit — > xp^ in (13.81) indicates that the 
above diagrams are evaluated in the approximation where each of the four-momenta of 
the various gauge bosons is collinear to one of the momenta p M of the external legs and/or 
soft. Where relevant, also the contributions of hard regions are taken into account (see 
Sect. 13. 2p . The factorizable two-loop terms (13. 8p are computed in Sect. 15.21 

The remaining two- loop NLL mass singularities, which we call non-factorizable (NF), 
are obtained by subtracting from the diagrams of type (13 .7ft the factorizable terms (13.81) 
and external self-energy diagrams. There are four sets of non-factorizable contributions 
that are associated with the first four diagrams in (13.71) . whereas the last diagram in (13.71) 
does not produce non-factorizable NLL terms. This is due to the fact that all NLL mass 
singularities resulting from the last diagram in (13 .7p originate from the region where the 
momenta of the gauge bosons Vi, V2 are soft and are thus included in the factorizable part 
of this diagram. 
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The non-factorizable terms associated with the first diagram in (13.71) read 



M 



NF,A 



EE E 

i=X 3=1 v m =A,Z,W ± 



-E 







(3.9) 



The non-factorizable contributions associated with the second diagram in (13.71) are 



M 



NF,_B 



EE E 







(3.10) 



*2 "''i 



The non-factorizable contributions associated with the third diagram in (13 .7p are 



M 



NF,C 



inn n 

EEE E 

z i=l 2=i fc=i y m =A,Z,W± 








(3.11) 



<-0 



Here the first subtraction term represents an external-leg self-energy insertion and the 
symmetry factor 1/2 is introduced in order to avoid double counting of the terms resulting 
from the permutation of the external legs j and k. The non-factorizable contributions 
associated with the fourth diagram in (13.71) read 



M 



NF,D 



EE E 

i=l j=i v m =A,Z,W ± 

3Ti 





FT TV 



14 



As indicated by the limits q% — > xpf, q± — > in (13.9p ~ fl3.12p . we extract the NLL mass 
singularities that appear in the above combinations of diagrams in the regions where the 
momentum of the gauge boson V\ is soft and the momentum q% of the gauge boson V2 
is collinear to the external momentum pi*. The terms (j3.9D — (j3.12j) are free of singularities 
in the soft limit x — > since all soft-soft singularities are included in the factorizable 
parts that are subtracted. 

The subtraction of the factorizable terms must be performed without double-counting 
of diagrams and regions, i.e. the terms that are subtracted in f l3.9p - fl3.12p must correspond 
exactly to the ones that are included in (13. 8p . This correspondence is not obvious at 
first sight, since certain diagrams appear a different number of times or with different 
symmetry factors in (JES) and (l3~9l)-( riT2jl . This is due to the fact that in (l3T9l)-( l3~T2|) the 
contributions of certain diagrams are split into various terms that result from different 
singular regions, whereas in (13.81) every diagram includes the contributions of all singular 
regions. For instance, let us consider the last diagram in (13. 9p . This diagram appears 
only once in (13. 8p but is subtracted four times in (I3.9I) - (I3.12I) . These four subtraction 
terms can be rewritten as 



where the last term in (13.131) corresponds to the two terms that are multiplied with a 
symmetry factor 1/2 in (13.111) . As can easily be seen in (I3.13p . these three contributions 
are associated with the three non-overlapping singular regions of the diagram that give 
rise to NLL terms0 and their sum corresponds to the complete contribution included in 
(13. 8p . Similarly, one can verify that all other subtraction terms in ( 13.9p -( !3~T2l) correspond 
exactly to the factorizable contributions included in ( 13.81) . 

In Sect. H]we will prove that the non-factorizable contributions (I3.9j) ( T37T21) cancel. 

3.2 Soft— collinear approximation for virtual gauge-boson exchange 

As discussed in the previous section, mass singularities appear when the external 
fermions emit virtual gauge bosons in the soft and collinear regions. Here we introduce a 

9 The region qg — > xp^f , q£ — > does not give rise to NLL mass singularities. 





(3.13) 
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soft-collinear approximation that describes the gauge-boson-fermion-antifermion vertices 
in these regions. This approximation permits to derive the soft and collinear singularities 
of generic n-fermion amplitudes in a process-independent way. 

Let us first consider an n-fermion diagram involving the exchange of a virtual gauge 
boson V = A, Z, W between the external leg i and some other (external or internal) 
lineH 

* = E G^(-q, Pl + g)^^ie 7 %V(^ «*)■ ( 3 - 14 ) 
^ (Pi + q) 2 1 

Yjp' 

Here G^,~'(—q,Pi + q) represents the truncated Green function corresponding to the in- 
ternal part of the diagram, which is depicted as a grey blob. The contraction with the 
spinors of the external fermions j = 1, . . . , n with j ' ^ i is implicitly understood. The loop 
momentum of the gauge boson is denoted as q and the propagator that connects the gauge 
boson lines and V 11 has been omitted. Commuting the Dirac matrices associated with 
the fermion propagator and the gauge-boson-fermion vertex and using the massless Dirac 
equation we have 

tyi + rf)iey«(Pi, Ki) = -e [2(pi + q)» - 7^] u{Pi, (3.15) 

In the collinear limit q^ — *> xp% the $u(pi, /%) term vanishes owing to the Dirac equation. 
Thus, we obtain 




lim^ m 5 i = ffi(-g, Pi + q)u( Ph 2 f ^ + qY , (.3.16) 

i^^Pi y ' (Pi + q) 

where we have introduced the shorthand G~(—q,Pi + q)lj = Y^tf. G~^^*(—q, P i + q)I^ Vi - 
The approximation (13.161) is applicable also to the case where the gauge boson V becomes 
collinear to another leg j ^ i. In this case, the term 7^ on the right-hand side of (13.151) 
is contracted with another soft-collinear factor (pj — q)^ resulting from the coupling of 
the gauge boson V to the leg j and 

lim ( Pj - q)jfi = (1 - x)xf j = 0. (3.17) 

Similarly, (I3.16P applies also to the case where the gauge boson V splits into two gauge 
bosons, V and V", with V being collinear to an external leg j 7^ i and V" soft. Here, 
the combination of the soft-collinear factor associated with the external leg j, the triple 
gauge-boson vertex and the term 7^ on the right-hand side of (I3.15P yields 

lim (pj - qY {-Ig^q^n + + 9»n"%'\ 7^ = °- ( 3 - 18 ) 



10 For the moment, we consider the case where the external particle i is a fermion. The generalization 
to antifermions is discussed below. 
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For the multiple emission of collinear gauge bosons V\ . . . V n by an external fermion i we 
obtain 



V£ n . . . Vf 1 v£ n vf 1 

(Pt + <?«) 2 (p» + <?i) 2 

(3.19) 

where = qi + . . . + qj. This approximation is also applicable to all cases where one of 
the gauge bosons V\, . . . , V n becomes collinear to one of the other external legs j ^ i and 
all remaining gauge bosons are soft. Therefore (13.191) provides a correct description of the 
gauge-boson-fermion couplings in all regions that are relevant for our NLL analysis. 

Also for the case where the external line % represents an antifermion, the emission of 
collinear gauge bosons produces factors —2el^ h (pi + qkY k and, apart from the obvious 
replacement of the spinor u(j>i, k^) by v(j)i, Kj), we obtain exactly the same formula as in 

fl3H). 

The soft-collinear approximation (13.191) permits to replace the Dirac matrices associ- 
ated with each gauge-boson emission by simple four- vector factors —2(pi + q~k) flk - In the 
soft limit, qj — > 0, these factors correspond to the well-known factors ~2p^ k that are used 
to derive soft singularities in the eikonal approximation. The soft-collinear approxima- 
tion (I3.19P can be regarded as an extension of the eikonal approximation that permits to 
describe the emission of gauge bosons in the soft and the collinear regions. When applied 
to the one- and two-loop factorizable contributions (I3.5P and (13. 8p . this approximation 
permits to factorize the mass singularities from the n-fermion Born amplitude explicitly. 

We note that the soft-collinear approximation is not applicable in the case where NLL 
two-loop contributions arise as a combination of logarithmic singularities originating from 
soft-collinear and UV regions. In particular, for the two-loop factorizable terms of the 
type 



(3.20) 



V' 



where a soft-collinear singularity resulting from the exchange of the gauge bosons V (and 
V) appears in combination with an UV singularity resulting from hard particles in one- 
loop subdiagrams, the approximation (I3.19P can be applied only to the vertices that occur 
outside the UV-divergent one-loop subdiagrams whereas for the vertices and propagators 
inside the one-loop subdiagrams we have to apply the usual Feynman rules. 

In the case of the last two diagrams in (13.201) this approximation is not sufficient to 
eliminate the chain of Dirac matrices along the external fermionic leg i. However, this 
chain of Dirac matrices can be contracted with the spinor of the external fermion by 
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means of a simple projector. Let us illustrate how this is achieved for the second diagram 
in([S2ni), 




-2e/J„ .G8$ ( Pl )X v ^ (p h Pj )u(pi, m) 



(3.21) 



with 



X v ^(p hPj ) 



(2n) D (g2_ M 2 )(pi + g)2fe ._ g) 2 • 



(3.22) 



Here the factor (pj — q)^ comes from the soft-collinear vertex along the leg j, and iT M 
represents the one-loop vertex. The truncated Green function G-i(pt) in (13.211) represents 
the internal part of the diagram, where the momentum q of the gauge boson V is set to 
zero according to the definition of the factorizable part. The contraction with the spinors 
of the external fermions j — 1, . . . , n with j '• ^ i is implicitly understood. If one contracts 
this Green function with the spinor of the fermion i one simply obtains the on-shell Born 
matrix element, 

K i )=M?' M "«' i "* n . (3-23) 

The chain of Dirac matrices occurring in ( 13.221) can be projected on the Dirac spinor 
using0 



X VVl ^(p i ,p j )u(pi, Ki) 

so that one obtains 

G^X^ip^p^uipi,^) 



Tr 



X V(pi ' P 'i(pi,Pj)u) Ki f'ip'j u(pi,Ki) 



(3.24) 



2piPj 



Tr 



The same trace projector can be applied to the last diagram in (13.201) . 



(3.25) 




2e 2 I^J^ j G^p l )X^(p l ,pM^^, 



(3.26) 



11 This identity can easily be verified by means of the general decomposition 

N 

fe=0ii=j,j hk=i,j 

the Dirac equation, rfiufjn, ni) — 0, the anticommutation relation = ZpiPj, the identities yf\ = 

tf? = 0, and the fact that, for massless fermions, the subamplitude X Vvi ^i(pi,pj) contains only Dirac 
chains with even numbers 2k of Dirac matrices. 
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where 



A . D f d D q (p J -qr^, l + i)i^(p l + q)^ i + 

MD J (27T)V {q 2 _ M 2 } [{pi + q)2] 2 {p . _ q) 2 



(3.27) 



and iY^^i represents the fermionic self-energy. Again we can project (I3.27P on the Dirac 
spinor using 



3.3 Ultraviolet singularities 



1 



Tr 



x^{p h p 3 )u) Ki i>4 3 



Mr 



(3.28) 



Let us first discuss our treatment of UV singularities at one loop. The UV singular- 
ities appearing in the bare loop amplitudes are cancelled by corresponding singularities 
provided by the counterterms. These cancellations give rise to logarithmic contributions 
of the form 

1 / Q 2 V 1 (Q 2 \ 




M 



oop , 



111 



Mr 



(3.29) 



oop / 



where the first and the second term between the brackets result from bare loop diagrams 
and counterterms, respectively. Here is the scale of dimensional regularization and 
we have factorized the term (pv/Q 2 ) e that we always absorb in a e [see (12.71) ]. The re- 
normalization scale is denoted as /ir, and pi oop represents the characteristic scale of the 
UV-singular loop diagram. This latter is related to the momenta of the lines that enter 
the loop. 

At one loop, since in the high-energy limit all combinations of external momenta 
are hard, the characteristic scale of UV-divergent diagrams that contribute to truncated 
n-fermion Green functions is always of the order p 2 
all logarithms that result from UV singularities 
counterterms. To this end, in our calculation we perform a minimal subtraction of all UV 
poles that appear in the bare loop diagrams and in the counterterms. The combination 
of these subtracted terms corresponds to 



,„, p ^ Q 2 . This permits us to isolate 
i.e. terms of the type (13.291) . in the 




Q 2 



l 



oop , 




(3.30) 



and is obviously equivalent to (13. 29ft . As a result of the minimal subtraction the loga- 



rithmic contributions e 1 [{Q 2 /p 
characterized by a hard scale pi 



loop/ 



oop 



— 1] originating from bare loop diagrams that are 
Q 2 vanish. 

Thus, at one loop we can restrict ourselves to the calculation of the mass-singular 
bare diagrams and the counterterms. The minimal subtraction of the UV singularities 
appearing in these two types of contributions permits us to ignore any other bare diagram 
that produces UV singularities. 

At two loops, pure UV singularities produce only NNLL contributions since every 
UV-singular loop produces only a single- logarithmic term. The only NLL two-loop terms 
resulting from UV singularities are combinations of one-loop UV logarithms with one- 
loop double logarithms resulting from soft-collinear gauge bosons. These terms originate 
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from one- loop insertions in the one-loop diagrams (13. 2p . Here, as in the one-loop case, we 
perform a minimal subtraction of the UV singularity, such that the logarithms associated 
with one-loop sub divergences from hard subdiagrams with /xf oop ~ Q 2 are completely 
isolated in the counterterms. As a consequence, the UV contributions associated with 
one- loop insertions in the internal part of the one- loop diagram (13.21) become irrelevant. In 
particular the UV logarithmic contributions resulting from the two-loop ladder diagrams 
of type 




(3.31) 



in the region where V\ is soft and V2 is hard vanish. Minimal subtraction of the UV 
singularities permits us to use the soft-collinear approximation for all gauge bosons Vi, V 2 
in the above diagrams despite of the fact that this approximation is not appropriate 
to describe the hard region and can in principle produce fake logarithms of UV origin. 
Instead, the diagrams of the type (I3.20p . which result from the insertion of one- loop 
subdiagrams in the lines that are not hard ( P f oop <C Q 2 ) in (I3.2p . give rise to non-negligible 
NLL contributions of UV type. These UV contributions are correctly taken into account 
in our calculation as explained in Sect. 13.21 



4 Collinear Ward identities and cancellation of non-factorizable terms 

In this section we prove that the non-factorizable subset of mass singularities, i.e. the 
one-loop terms (13. 6p and the two-loop terms (I3.9p - (l3.12p . vanish. The proof is based on 
the collinear Ward identities that have been derived in Ref. [7]. 

Let us start with the non-factorizable one-loop terms (13.61) . This contribution can be 
written as 

m7 = ± E E & d I Zl 77, J, ,„ ^ ft* , 2( Pt - q y 

i=l V=A,Z,W± ip'. 



{2-K) D {q 2 -M^){p i -qf q ^„ : 

x { G^\ q , Pi - q )u( Pi , *) + EE ^^^^^ iW }^«. 

(4.1) 



... <p) 



where we have factorized the two propagators that appear in all three diagrams in (13.61) . 
Here u(pi, Ki) is the spinor of the external fermion i and we have introduced the abbrevi- 
ation 



V u r V u 



for the tree subdiagrams that are associated with the first two terms on the right-hand 
side of (I3.6p . In (14. 2 j) the contraction with the spinors of the external fermions j — 1, . . . , n 
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with j ' ^ % is implicitly understood, and the incoming momenta associated with the lines 
V/j, and i are q and Pi — q, respectively. 

The combination of tree subdiagrams (14. 2 p fulfils the collinear Ward identities [7] 



lim q^G^~- li (q, p { - q)u(p h Ki) = ^M 

P — >xn 



tpi—ifi—tpn T V 
61 <fi'M- 



(4.3) 



Using the charge- conservation relation (I2.30P and 

lim 2 ^ + ^ = l 
we can rewrite the identities (14.31) as 



for j ^ i, 



q^^xp 



(4.4) 



0, (4.5) 



- fi 



where the expression between the curly brackets is identical to the one that appears on 
the right-hand side of (14.11) . These collinear Ward identities can be represented diagram- 
matically as 

JCf— < -E 1 



lim x 




3 = 1 



0. 



(4.6) 



where the contraction with all fermionic spinors, including u(pi, Ki), is implicitly under- 
stood. 

The cancellation of the non-factorizable terms 1 14. ip is simply due to the fact that in the 
collinear limit g M — > xp^ the four- vector (p t — q)^ on the right-hand side of (14. Xp becomes 
proportional to the gauge-boson momentum g M and its contraction with the expression 
between the curly brackets vanishes as a result of the collinear Ward identities (14.51) . 

Let us now consider the two- loop non-factorizable terms (13.9 j) — (13. X2p . The contribution 
( 13T9D yields 



M 



NF,A 



V"V" 2(4-D) /" d D gi f d D q 2 8e 3 (pi - q x )(pj + q x ) 

2^2^ 2^ 2^ 



*=1 i- 1 V m =A,Z,W± V ' ip 1 ' ip'. 



X 



(Pi - 9l) 2 (Pj + Ql) 2 {Pi ~qi- Q2) 2 q^Oq^xp? 
X - (g 2 ,Pi - ft - q2MPi, Ki) 



(2n)°J (2ir)° (qf - M Vi )(q : 
lim \im (p i -q l -q 2 y 



+ 



2(Pj + qi + Q2) fi x - ,^1...^...^...^^^ 



(Pj + 91+92) 



E^ 



fc =i (pk + 92) 



ipi...^ ...tp j ...ip k ...tp n y 



jV 2 jVi rVx 

V' k <Pk I ViVi f'jfj f'ifi 



0. (4.7) 
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This cancellation can easily be verified by means of the collinear Ward identities (14.51) 
by observing that in the soft-collinear limit q± — > 0, q 2 — > xp? the four- vector (pi — q\ — 
q 2 y L tends to (1/x — l)q 2 an d the expression within the curly brackets in (14. 7B becomes 
equivalent to the one in (14.51) . Similarly, for the contributions (I3.10P and (13.111) we obtain 



M 



NF,B 



2^ 2^ 2^ z_/ 



i=l 3=1 V m =A,Z,W± ip> ip"ip'. 

J7=2 lr lr J 



frr)°J (2tt)^ (q 2 - M^){ql - M 2 ) 



ii m „ lim >< - ^y 



(Pi - Q2) 2 (Pj + qi) 2 (Pi -qi- Q2) 2 tf-o^-asp? 
x <^ Gji * [q 2 , Pi - qi - q 2 )u(pi, Ki) 



, 2(pj + ft + g 2 )„ ^ Vl ... v ^... V '!... Vn v 2 

+ fe + gi + g 2 ) 2 ^ 



and 



1 n n n r AD „ r AD „ i 

1 \ - \ - \ - \ - \ - 2(4-0) f a qi f a q 2 J 



(2^7 (2tt)° (qf - M 2 J 



x 



8e 3 (p k - gi)(pj + q x ) 



i im >< - 



(ql - M 2 2 )( Pi - q 2 ) 2 ( Pj + gi ) 2 (p k - qi f rf-orf-*^ 



x \ G [ J 2 -\q 2 ,Pi - q 2 )u(pi, k { ) 



+ 



2(Pj + gi + gg) 
(Pi + 9i + 92) 



pT v 2 



+ 



2(pk - gi + g 2 ) A 
(Pik - ?i + g 2 ) 2 



+ E 



tt (Pi + Q2) 2 



eI V2 M Vl I n V 7 

f'lfl f <Pj<Pj ^'fcVfe f'ifi 



rV 2 



Finally, for the contribution (13 . 1 2[) we have 



0. 



(4.9) 



M 



NF,D 



EE E E/4 1 



2(4-D) 



d% /■ d w g 2 



4ie 2 5r 2 ^ 



V1V2V3 



(2*)° J ^(gf-M^-M^) 



x 



X 



i im J im ufe ~ ^) M2 (Pj - qi) 



(gf - M^)(pi - g 2 ) 2 fe - gi) 2 tf-otf-** 
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(Pj + Q2) 



(4.10) 



where q% = q\ + q%. In the soft-collinear limit we observe that 

£ m i im M ~ 93)^ (Pi ~ <?ir 9 win ~ + < 2 3 fe + - <'(g3 + 9i 



9i— *0 g£-«cp? 



U'2 



fj-.i 



(4.11) 



and again the contraction of this four- vector with the expression between the curly brack- 
ets in ( I4.10p cancels as a result of the collinear Ward identities (14.51) . 



5 Factorizable contributions 



In this section, we present explicit results for the one- and two-loop factorizable contri- 
butions defined in Sects. I3.1.2l and l3.1.4[ These are evaluated within the 't Hooft-Feynman 
gauge, where the masses of the Faddeev-Popov ghosts u A , u z , u w and would-be Gold- 
stone bosons x,^ read M u a = M A = 0, M x = M u z = M z , and M ± = M u w = M w . 
Using the soft-collinear approximation introduced in Sect. 13.21 we express the factorizable 
contributions resulting from individual diagrams as products of the n-fermion Born am- 
plitude with matrix-valued gauge couplings and loop integrals. The definitions of these 
loop integrals are provided in App. [A] The integrals are computed in NLL accuracy, 
and the result is expanded in e up to C(e 2 ) at one loop and O(e ) at two loops. The 
UV poles are eliminated by means of a minimal subtraction as explained in Sect. 13.31 
such that the presented results are UV finite. All loop integrals have been solved and 
cross-checked using two independent methods: an automatized algorithm based on the 
sector-decomposition technique [33] and the method of expansion by regions combined 
with Mellin-Barnes representations [34]. 

5.1 One-loop diagrams 

The one-loop factorizable contributions originate only from one type of diagram J^l 

= -M £ I^I^D (M Vl ; rij ). (5.1) 

V 1 =A,Z,W ± 

3 

The corresponding loop integral Do is defined in flA.61) and to NLL accuracy yields 
A)(M W ; r«) N ^ L —L 2 - \l\ - ~L 4 e 2 + 2 (2 - Z y ) (l + \h\ + ~L 3 e 2 ) , 

A)(M Z ; r y ) N = L A)(M W ; r y ) + h (2L + 2L 2 e + L 3 e 2 ) , 

£>o(0;r y ) = L -2e- 2 - 2 (2 - e~\ (5.2) 
12 The Z-loop diagrams depicted in this section are understood without factors (a e /47r) ! . 



23 



where the UV singularities 

C(Mw; r l3 ) ^ D™(M Z ; r y ) ^ D UV (0; r y ) N ^ L 4c" 1 (5.3) 
have been subtracted. 

5.2 Two-loop diagrams 

The two-loop NLL factorizable terms (13.81) involve fourteen different types of diagrams. 
The diagrams 1-3, 12, and 14 in this section give rise to LLs and NLLs, whereas all other 
diagrams yield only NLLs. The loop integrals associated with the various diagrams are 
denoted with symbols of the type Dh(rni, . . . , m?; Tm, ■ ■ .) and depend on various kine- 
matical invariants and masses m$. The symbols are always used to denote generic 
mass parameters, which can assume the values rrii = Mw, M%, m t , or = 0. Instead 
we use the symbols Mj to denote non-zero masses, i.e. M, = M\y, Mz, m t , Mh. The inte- 
grals are often singular when certain mass parameters tend to zero, and the cases where 
such parameters are zero or non-zero need to be treated separately. For every integral 
we first evaluate D/j(Mw, . . . , Mw; r^i, ■ ■ .), i.e. the case where all mass parameters are 
equal to Mw- The dependence of the integral on the various masses is then described by 
subtracted functions of the type 

AD h (mi, rrij] r M , . . .) = D h (mi, . . . , mj] r M , . . .) - D h (M w , . . . , M w ; r M , ...). (5.4) 
Diagram 1 



MY j = KfXJM" 1 =M ° £ ijHYHjHfD^My^My^l (5.5) 

Vi,V 2 =A,Z,W ± 

~3 




where the loop integral D\ is defined in ( 1A.6|) and yields 



D 1 (M W , M w ; r y ) = L \L A - H (2 - l tJ ) L 3 , 
AD l (M 1 ,M 2 ;r ij ) N ± L -^hL 3 , 



AD 1 (0, M 2 ; rij ) N i^ L 2L 2 e- 2 + ^LV 1 + ^L 4 - (2 - /<,-) (ALt- 2 + AL 2 e~ l + 2L ; 



-l 2 {ALe' 2 + 8L 2 e~ 1 + 8L 



AD 1 (M 1 ,0;r ij )"= —hL 3 , 

ADi(0, 0; r y ) = L e~ 4 - ^L 4 + (2 - ^e^ 3 + ^L 3 ) . (5.6) 
Here the UV singularities 

<(M 1)m2 ;r,)1 L -4LV 1 -^ 1 
D^ v (0,m 2 ; rij ) N ^ L -8e" 3 (5.7) 

have been subtracted. 
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Diagram 2 




M J2 lf 2 lf 1 lJ 1 lJ 2 D 2 (M Vl ,M V2 ;r 

Vi,V 2 =A,Z,W ± 



(5i 



where the loop integral D 2 is defined in (IA.6I) . This integral is free of UV singularities 
and yields 



D 2 {M w ,M w ; ri 
AD 2 {M 1 ,M 2 ;r i 
AD 2 (0,M 2 ;r, 



NLL 



NLL 



NLL 



•L 4 — — (2 — lij) L 3 , 

~L\- 1 - 7 -L 4 + (2 - kj + h) ( 2L 2 e- 1 + ^-L 



AD 2 {M X , 0; r l3 ) W -h 3 e^ - + (2 - l {j + h) (2LV 1 + ^L 3 



AD 2 (0, 0; nj ) ^ L 6~ 4 - + (2 - iy) (V 3 + ^L 3 ) 



Diagram 3 



(5.9) 




-i^M e v ^ v U^I^I^D^M Vll M V2 ,M V3 -r l0 ), (5.10) 



where the e-tensor is defined in (12.181) . The loop integral D 3 is defined in (I A. 6 1) and yields 

D 3 (M w ,M w ,M w;rij ] 



NLL l L 4_/ 3 



6 



^3 



AL>3(M 1; M 2 , M 3 ; r y ) N = L -~ ft + Z 3 ) L 3 , 



AD 3 (0,M 2 ,M 3 ;n 



NLL 



^L 4 + (2 - Zy + Z 3 ) (L^ + h 3 ). 



AD 3 (Mi, 0, M 3 ; r tJ ) N ± L -i ft + Z 3 ) L 3 , 

AL> 3 (M 1; M 2 , 0; r y ) = L -^V 1 - ^ - QLe~ 2 - (2 + Z y - Z x ) L 2 e _1 

/11_5Z, 5ZA 
3 3 3 / ' 



(5.11) 



where the UV singularities 

DY v (mi, m 2 , M 3 ; r„ 



-3LV 1 - 2L 3 , 



(5.12) 
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have been subtracted. 



Diagram 4 




* =-m e n 2 ipirij i D 4 {M Vll M V2 

Vi,V 2 =A,Z,W ± 

-j 

where the loop integral _D 4 is defined in ( 1A.6|) and yields 



D 4 (M w ,M w ;r, 



N = L -L 3 



AD 4 (Mi, M 2 ; Tij) N ^ L 0, 



AD 4 (0, M 2 ;r, 
A J D 4 (M 1 ,0;r JJ ) N i± L 0, 



AD 4 (0,0;r 



Here the UV singularities 



NLL _ e _ 3 _ 1^3 



NLL r 2 -1 



NLL „ _ 3 



have been subtracted. 



Diagram 5 




-M E iy'lM'lf'DsiMv^My, 

Vi,V 2 =A,Z,W ± 



where the loop integral D 5 is defined in ( 1A.6|) and yields 

D 5 (M w ,M w ; rij ) "=-h 3 , 



AD 5 (M 1 ,M 2 ;r ij ) N ± L 
AL> 5 (0, M 2 ;r ; 



NLL 



-2Le- z - 2LV 1 - L a , 



AL> 5 (M 1 ,0;n i ) N ^ L 0, 



AD B (0,0;r. 



W e~ 3 + -L 3 . 



Here the UV singularities 



nUV/ Ti/r \ NLL r 2 -1 ^ 

L> 5 {M 1 ,m 2 ;r ij ) = —L e 



-£ 3 , 
3 ' 



£> 5 uv (0,m 2 ;r. 



have been subtracted. 
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Diagrams 6 



M 6 / j = 




+ 




102 



2^ 



v 1 y 2 y 3 y i = 



3,Va=A,Z,W± 

where the loop integral D e is defined in ( 1A.6|) and yields 

20 

As(M w , M w , M w , M w ; r -) N = L y 



k(M Vl ,M V2 ,M V3 ,M V4 ; 

r ij)i 

(5.19) 



AD 6 (M 1 ,M 2 ,M 3 ,M 4 ;r lj ) N = 0, 



AD 6 (0, M 2 ,M 3 ,M l 



A/AiLU,.!)..!/,..!/;:, 

A J D 6 (M 1 ,M 2 ,0,M 4 ;r u ) N = L 



NLL M| + Mf 

" 2M| 



-16Le- 2 



. 8L 2 e -l + 1^ L 3 

3 



\ NLL n 

4;^) = 0, 



AD 6 (M 1 ,M 2 ,M 3 ,0;r l ,) N i± L ^|±M [_i 6Le - 2 



2Mf 

AD 6 (0, M 2 , M 3 , 0; r y ) N ^ L yLe~ 2 + y LV 1 - ^ 
Here the UV singulariti* 



3LV + 15 



10 

+ T 



i 3 



ties 



ClM.m.ma,!^,) N ^ L ^LV 1 + ^L 3 , 
D 6 UV (0, ma, m 3 , M 4 ; r «) N ^ L y LV 1 + y-^ 3 

TTln + 771a 

J £ ± —I 

2Mi [ 

D^(M 1} m 2 ,m 3 ,0;r l3 ) ^ y^V 1 + y-L 3 



-lQe- 3 + 8L 2 e- 1 + -, 
3 



+ 2 7,9 -16e- 3 + 8LV 1 + — L 3 



D 6 uv (0,m 2 ,m 3 ,0; 



\ NLL 20 



2M\ 



(5.20) 



- . „• (5-21) 

have been subtracted. We observe that the loop integrals associated with A—Z mixing- 
energy subdiagrams give rise to the contributions 

AAs(0, M W ,M V 



/If 2 

fw,M z ;r y ) = ^ 



which depend linearly on the ratio M^/M^. 



-16Le- 2 - 8LV 1 + ^L 3 



(5.22) 
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Diagram 7 




4m 



E iM 3 E 

VuV 2 ,V 3 =A,Z,W ± V=A,Z,W ± 



£ V 1 V2V £ V3V 2 V 



x{D- l)D 7 (M Vl ,M V2 ,M V3 ;r t3 ), 
where D = 4 — 2e. The loop integral D 7 is defined in ( 1A.6|) and yields 



D 7 (M W , M w , M w ; r, 
AD 7 (M ll M 2 ,M 3 ;r i 



NLL 



NLL 



o, 

M 2 

AD 7 (0,M 2 ,M 3 ;r ii ) r m L _| 
AD 7 (M 1 ,0,M 3 ;r 4J ) N i± L 0, 



AD 7 (M 1 ,M 2 ,0;r ij ) N ^ L ^ 



^ M 2 
lj) ~ M 2 



-2Le~ 2 - LV 1 + -L 3 
3 



-2Le- 2 - LV 1 + -L 3 
3 



AD 7 (0,M 2 ,0;r y ) N = L 0, 



where the UV singularities 



D 7 uv (0,M 2 ,M 3 ;r,) N ^S 



D^(M 1 ,M 2 ,0;r ij ] 
have been subtracted. 



M| 
Mf 

NLL Mf 

M? 



-2e 



~ 3 + L 2 e _1 + ^L 3 



2 

— j 

3 

2 



-2e~ 3 + LV 1 + -L 3 
3 



(5.23) 



(5.24) 



(5.25) 



Diagram 8 



mt 3 




-e 2 v 2 M E ^M* E 



{l v * , I* } <s>iH D 8 (M Vl , M^ 2 , M V3 , M Vi ; r y ) , 



(5.26) 



where the curly brackets denote anticommutators and v is the vacuum expectation value. 
The loop integral D s is defined in ( 1A.6|) and yields 



M W D 8 (M W , M w , M w , M w ; r^) N = L 
U 2 x AIK(U...\l>..\I A .U l :r, 



NLL 



AL^M^M^M^ry) 



NLL 



M| 



-2Le~ 2 - L 2 e _1 
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M^AD 8 {M u M 2 ,0,M 4 ;n 



NLL 



0. 



A J D 8 (M 1 ,M 2 ,M 3 ,0;r li ) N ^ L ^ 



-2Le- 2 - LV 1 + -L 3 
3 



M^AD s (0,M 2 ,M 3 ,0;r, 



NLL 



'.I i 



0. 



Here the UV singularities 



D 8 uv (0,M 2 ,M 3 ,M 4 ;r, 



NLL 



Mf 

,1'V, w w w \ NLL 1 



D 8 uv (M 1 ,M 2) M 3 ,0;r ii 



M 2 



-2e~ 3 + LV 1 + -L 6 
o 

-2e~ 3 + LV 1 + -L 3 
3 



(5.27) 



(5.2* 



have been subtracted. The above diagram represents the only contribution involving 
couplings proportional to v, which originate from spontaneous symmetry breaking. 

Diagram 9 



M 9 / j = 




Vi,V A =A,Z,W± i2! i3 =H,x,^ ± 



X 



J 

D 9 {M Vl ,M* i2 ,M® i3 ,M V4 ; 



where the loop integral D 9 is defined in ( 1A.6j) and yields 

A,(M w ,M w ,M w ,M w ;ry) N ^ L ^L 3 , 



(5.29) 



9 

AD 9 {M 1 ,M 2 ,M 3 ,M 4 ;r ij ) N = 0, 
AD 9 (0,M 2 ,M 3 ,M 4 ;r,)^ L ^|±^ 



4Le" 2 + 2LV 1 - -L 3 
3 . 

ALe' 2 + 2LV 1 - -L 3 



AD B (0, M 2 , M a , 0; r y ) N ^ L ^e~ 2 + IlV 1 - ^L 3 . 



Here the UV singularities 

J D 9 uv (M 1 ,M 2 ,M 3 ,M 4 ;r i ,) N ^ L l 



(5.30) 



IT (0, M 2 , M 3 , M 4 ; r,,) ^ ^V* + ^L 3 + ^±^1 
D™( Ml , M 2 , M 8) 0; r y ) W itfe" 1 + ^L 3 + M±M 

o 



4e" 3 - 2LV 1 - -L 3 
3 



D 9 uv (0,M 2 ,M 3 ,0;r^) N ^^ 



4e~ 3 - 2LV 1 - -L s 
3 



(5.31) 



have been subtracted 
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Diagram 10 




-> E ^M 3 E {i v \i Vi } 

Vi,V 3 =A,Z,W ± *i 2 =/f,x,0 ± 



x D 10 (M Vl ,M*.,M V3 ; 



(5.32) 



where 



10 



£> 7 . 



(5.33) 



Also this diagram, which yields NLL contributions only through A—Z mixing-energy sub- 
diagrams, gives rise to a correction proportional to M^/M| originating from 
ADj(Q,M^±,Mz;rij). This correction cancels the contribution proportional to M|y/M| 
that originates from diagram 9. 

Diagram 11 

For the diagrams involving fermionic self-energy subdiagrams we consider the contri- 
butions of a generic fermionic doublet ^ with components = u, d. The sum over the 
three generations of leptons and quarks is denoted by and colour factors are implic- 
itly understood. Assuming that all down-type fermions are massless, rrid = 0, and that 
the masses of up-type fermions are m u = or m t , we have 



Ml 1 ''' 




V 1 ,V i =A,Z,W ± 



X 



El E E D nj o(M Vl ,mi 2 ,m i3 ,M v ^rij) 

' * * ' * fi to 19 to 



(5.34) 



where Du )m = — 4D 8 represents the contribution associated to the m^-terms in the nu- 
merator of the up-type fermion propagators, whereas the integral -Dn,o, which is defined 
in (lA.6p . accounts for the remaining contributions. This latter integral yields 



D 11>0 (M w ,M w ,M w ,M w ;r. 
ADu )0 (Mi, m 2 , m 3 , M 4 ; r, 
AD nfi (0, m 2 ,m 3 ,M 4 ;r. 



n. 



vn 



AD Ufi (M 1 ,m 2 ,m 3 ,0;r. 



NLL o 3 
— —J-i , 



NLL 



nll ml + ml r 

2M| 
nll m\ + m\ 



VI) 



2Ml 



]Le~ 2 + 4LV 1 - -II 



ILe^ + ALh- 1 -^ 3 
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AD nfi (0, M 2 , M 8 , 0; r,,) N ± L ^ + ^V 1 - \l\ 



AD 11)0 (0, 0,0,0;^-; 
where the UV singularities 



NLL 



_ 2e -3 _ ° L 3 

9 ' 



D^ (Mi,m2,m3,M 4 



NLL ^ r 2 -1 



) ' ill 



9 



nuv/ n m , f \ nll 4 2 i 8 3 m 2 + m 3 
i^ 110 (U, m 2 , m 3 , M 4 , r^-j — -l, e + -L h — 



D li y o (Mi,m 2 ,m 3 ,0;r i 

(0,7712,7713,0;^ 



nll 4 r2 _! 8 3 m^ + m^ 



2 r 



^e" 1 + + 



2MI 



- 3 - 4LV 1 



" 3 - 4LV 1 



have been subtracted. As a consequence of 

AD 1 i t0 (0,m 2 ,m 3 ,M 4 ;r ij ) = — - ADn )m (0, m 2 , m 3 , M 4 ; ?-„), 
AZ) 11)0 (Mi, 7712,7713,0;^) = - — -AD llim (M 1 ,m 2 ,m 3 ,0;r ti ), 

all terms proportional to the fermion masses in A4 z ,t3 cancel. 
Diagram 12 




M Yl I? 2 lFl?I?D 12 (M Vl ,M V2 

V 1 ,V 2 =A,Z,W ± 



where the loop integral D 12 is defined in (IA.6I) and yields 



D 12 (M W , M w ; r lk ) N ± L -L 4 - 2 (2 - L 3 , 
AD 12 (M U M 2 ; r lk ) N = L (2 Z x + Z 2 ) L 3 , 
AL> 12 (0, M 2 ; r ifc ) N ^ L 2L 2 e~ 2 + 2LV 1 + ^L 4 - (2 - l lk ) ^Le" 2 + 2L 2 e 



2-1 



Z 2 ( 4Le~ 2 + 6LV 1 + 4^L 3 , . 



14 
~3~ 



10 



AD 12 (Mi, 0; r lk ) == - -LV 1 - ^L 4 + (2 - Z jfc ) ( 2LV 1 + 

+ Zi (2L 2 e- 4 + h 
AD 12 (0, 0; r 4fc ) N ^ L 2e~ 4 - + (2 - l lk ) (a^ 3 + 2L 3 ) . 
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Here the UV singularities 



°L 3 , 
3 ' 



— »e 



(5.40) 



have been subtracted. Note that, to NLL accuracy, the above diagram does not depend 
on Tij and Tj k . 

Diagram 13 



Mf' ljk 




E e™ v H?Uj*l¥D l3 (M Vl ,M V2 ,M, 



V3 j > T"ik j ^ jkji 



V 1} V2,V 3 =A,Z,W± 



(5.41) 



where the loop integral -D13 is defined in (I A . 6 1) . This integral is free of UV singularities 
and yields 



D 13 (M W , M w , M w ; nj, r ik . 


r ik) 


NLL 


0, 






AD 13 (M 1 ,M 2 ,M 3 ;r ij ,r ik , 


r ik) 


NLL 


0. 






AD 13 (0,M 2 ,M 3 ;r ij ,r ik , 


r ik) 


NLL 


{hj 


— hk) 




AD 13 (M 1 ,0,M 3 ;r ij ,r ik , 


r ik) 


NLL 


(Ijk 


— hj) 




AD 13 (M 1 ,M 2 ,0;r ij ,r ik , 


r ik) 


NLL 


(hk 


- Ijk) 





3 
5 



(5.42) 



Diagram 14 



Ml 4 ' ijkl 




M E I^I^I^I^D 14 (M Vl ,M V2 ;r ij ,r kl ) 

Vi,V 2 =A,Z,W ± 



(5.43) 



where the loop integral D u is simply given by the product of one-loop integrals (15. 2p . 

D u (M Vl ,M V2 ;nj, r kl ) = D (M Vl ; r y ) D (M V;2 ;r kl ). (5.44) 
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6 Renormalization 



In this section we discuss the NLL counterterm contributions that result from the 
renormalization of the gauge-boson massesf^l 

M Vj0 = M y + f;(^y<5M<?, (6.1) 

and the electroweak couplings, 

m = to + £ (t) 1 89?, e = e + £ (^Se®, (6.2) 



1=1 7 1=1 7 

as well as from the renormalization constants associated with the wave functions of the 
external fermions k — 1, . . . , n, 



00 7 a f ^ 



The renormalized one- and two-loop amplitudes are presented in Sect. [7J 

6.1 One-loop contributions 

At one loop, the mass counterterms SMy^ can be neglected, since the gauge-boson 
mass terms in the Born amplitude give only contributions of order M v /Q 2 , which are 
suppressed in the high-energy limit, and the same holds for the contributions resulting 
from their renormalization. 

The electroweak couplings are renormalized in the MS scheme with an additional 
subtraction of the UV singularities as explained in Sect. 13.31 Assuming that the renor- 
malization scalJ^l /zr is of the order of or larger than Mw, this yields the counterterms 



5 (1) NLL _9il (1) 

2 e 




, (l) NLL _£ 1,(1) 

- 2e e 




(6.4) 



and the one-loop /i-function coefficients bi , b^ , and are defined in App. O The 
dependence of the counterterms (16.41) on the factor (Q 2 //xr) £ is due to the normalization 
of the expansion parameter a e in (16.21) . As explained in Sect. 13. 3[ the expressions (16.41) 
are obtained by subtracting the UV poles from the usual MS counterterms. Since the 
same subtraction is performed in the bare loop diagrams, the resulting renormalized 
amplitudes correspond to the usual MS renormalized amplitudes. The renormalization 
of the mixing parameters c w and s w can be determined from the renormalization of the 
coupling constants via (12.221) . In NLL approximation, this prescription is equivalent to 
using the on-shell renormalization condition, i.e. relation (12.211) . for the weak mixing angle. 

Here and in the following we assume that the Born amplitude Ai is expressed in 
terms of coupling constants renormalized at the scale //r = Q. As a consequence, the 
contribution of the counterterms (16.41) to the one- loop amplitude M. 1 vanishes. 



13 Note that in (|6.ip - (|6.3p we use the expansion parameter a c defined in (|2.7[) . 

14 We do not identify the renormalization scale /zr and the scale of dimensional regularization /id- 
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The only one-loop counterterm contribution arises from the on-shell wave-function 
renormalization constants SZ^ for the massless fermionic external legs. These receive 
contributions only from massive weak bosons, whereas the photonic contribution vanishes 
owing to a cancellation between UV and mass singularities within dimensional regular- 
ization. After subtraction of the UV poles we find 



SZ, 



(1) NLL 



E Ik Ik 

v=zw± 



- 1 



- I A I A 



(6.5) 



Finally, the one-loop counterterm for a process with n external massless fermions in NLL 
approximation is obtained as 



n -r 

WF _ a a \ " 1 f 7 (l) 



k=l 



(6.6) 



6.2 Two-loop contributions 

At two loops, the mass renormalization leads to non- suppressed logarithmic terms 
only through the insertion of the one-loop counterterms 5My^ in the one-loop logarithmic 
corrections. However, these contributions are of NNLL order and can thus be neglected 
in NLL approximation [28]. In this approximation also the purely two-loop counterterms 
that are associated with the renormalization of the external-fermion wave functions and 
the couplings, i.e. dzj^\ 5g^ and 5e^ 2 \ do not contribute. 

The only NLL two-loop counterterm contributions are those that result from the com- 
bination of the one- loop amplitude with the one- loop counterterms 5Zjp, 5g^ and 5e^\ 
The wave-function counterterms yield 



M 



WF 
2 



Ml 



F \ ~» 



k=i z 



(6.7) 



The wave-function renormalization constants bZ k and the unrenormalized one-loop am- 
plitude Mi are given in (I6.5P and in (ID. II) . respectively. In NLL approximation the 
one-loop counterterms Mi do not contribute, and only the LL part of the one-loop 
amplitude M\ is relevant for (16. 7p . This is easily obtained from (ID. 1[) by neglecting the 
NLL terms depending on and lz, and using global gauge invariance f!2.30|) as 



e 2 M\^ l -M,±{ £ e 2 /r/r^o(M w ;-Q 2 ) + e 2 /^A/J (0;-g 2 )^ (6.8) 



LL 1 



i=l I V=A,Z,W ± 
n 



2-^0 E 

z i=i 



D (M W ; -Q 2 ) + e 2 Q 2 AD (0; -Q 2 ) . (6.9) 



In the second line we used the identity (12.161) . Note that only the LL contributions of D Q 
and AD contribute in (I6.8P and (16.91) . 

The remaining NLL two-loop counterterms result from the insertion of the one-loop 
coupling-constant counterterms (16 .4p in the LL one-loop amplitude (16.91) and read 



e 2 Ml K N ± L -- 
2 2e 




- 1 



M,E< 



i=l 



Yi 



+ g 2 2 b?a 



D (M W ; -Q 2 ) 
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+ e 2 6WQ, 2 A J D (0;-g 2 )|. (6.10) 

The various one- loop /3-function coefficients in (16. 101) . op, and b£\ are denned in 
App.O 

7 Complete one- and two-loop results 

In this section we combine the unrenormalized and the counterterm contributions pre- 
sented in Sects. [5] and [6] and provide explicit NLL results for the renormalized one- and 
two-loop n-fermion amplitudes. As we have seen in Sects. and El the NLL corrections 
factorize, i.e. they can be expressed through correction factors that multiply the Born am- 
plitude. Moreover, we find that the two-loop correction factors can be entirely expressed 
in terms of one-loop quantities. In this section we concentrate on the results. Details on 
the combination of the various contributions can be found in Apps. 03] and [El 

7.1 Renormalized one-loop amplitude 

The renormalized one-loop matrix element for a process with n external massless 
fermions is given by 

M x = M\ + M™*, (7.1) 

where 

inn 

^= EEK (7.2) 

Z 1=1 3=1 

represents the bare one- loop contribution, which is given by the factorizable part ( 13.51) . 
and Mi is the wave-function-renormalization counterterm (16. 6p . Using the explicit 
results presented in Sect. 15.11 we can write (more details can be found in App. [Dj) 

Mi N = L M [Ff w + AF? m + AF-f] . (7.3) 

Here the corrections are split into a symmetric-electroweak (sew) part, 

inn 

^T CW = -^EE E if///(6,M w ;-r«), (7-4) 

z i=i j=l v=A,Z,W± 

which is obtained by setting the masses of all gauge bosons, A, Z and W , equal to Mw 
in the loop diagrams, an electromagnetic (em) part 

A*T = ~\ E E ft} 0; -r„), (7.5) 

Z t=l 3=1 

resulting from the mass gap between the W boson and the massless photon, and an 
M z -dependent part 

inn 

Aif = --££/f/;A/(e,M z ;-^), (7.6) 

Z 1 = 1 3=1 

3=£i 
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describing the effect that results from the difference between Mw and M%. For the 
functions 7, including contributions up to the order e 2 , we obtain 

7(e, M W ; -r«) N ^ L -7 2 - ^7 3 e - iiV + (3 - 2iy) (l + ~L 2 e + IlV) + 0(e 3 ), 

J(e, M z ; N = L 7(e, M w ; -r y ) + Z z (27 + 27 2 e + 7 3 e 2 ) + 0(e 3 ), 
7(e, 0; -n 3 ) ^ -2e- 2 - (3 - 2Z ii )e" 1 , (7.7) 
and the subtracted functions A7 are defined as 

A7(e, m; -r y ) = 7(e, m; -r^) - 7(e, M w ; -r y ). (7.8) 
In LL approximation we have 



7(e, m; -r i3 -) = 7> (m; r^), 



(7.9) 



such that we can replace 7? and A7) by 7 and A7, respectively, in (16. 8p . (16.91) . and 
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7.2 Renormalized two-loop amplitude 

The renormalized two-loop matrix element reads 



M 2 = M\ + Mf F + M P 2 R , 



(7.10) 



where 



II IL r -1 



= 1 3=1 



2 



I U n l i i n 

+ =; E ^ + E A*?** + --MP* + 5 E Ml 4 

1 m=6 k=l I b 8 1=1 



(7.11) 



represents the bare two-loop contribution, which is given by the factorizable terms (13.81) . 
and M^ F and M 2 R are the wave-function- and parameter-renormalization counterterms 
given in (16.71) and (16.101) . respectively. As shown in App. [Ej the two-loop amplitude (17.101) 
can be expressed in terms of the Born matrix element M and the one-loop correction 
factors fl73)- (T76]) as 

M 2 N = L M l i [Ff™] 2 + F?™AF* m + i [AFf 11 ] 2 + Ff^AF? + AFfAFf 11 



+ G s 9 ew + AG 



cm I 

2 r> 



(7.12) 



where the additional terms 



1 n 

e 2 G sew = _ ^ 
^ i=l 
1 " 

2 



AG: 



cm 



J(e,M w ,/4), 

E Q 2 ] ^ i A J ( e ' °> ^) - AJ ^ °> M w)] + & Qed A J(e, 0, M 
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w) 



(7.13) 



contain one-loop /3-function coefficients, defined in App. [UJ and the combinations 



J(e,m,/4) 



I(2e,m,Q 2 )- (%\ I(e,m, Q 2 



A J (e, m, /4) = J (e, m, /4) - J (e, M w , /4) 



(7.14) 



of one-loop /-functions (17. 7p and (17.81) for m = M w , M z , 0. The relevant J-functions read 
explicitly 



J(e,M w ,// R ) 
A J(e, 0, M 



2 ^ NLL 1 ^ 3 



3 

NLL 3 

= 2 C 



/ T 



- 3 + 2Le- 2 + LV 1 + 0(e) 



AJ(e,0,/4) - AJ(e,0,M, 



where 



NLL 



"MR 



-2e" 



e - 1 (/ MR -2L) + /, R L--/ 2 R 



/<R 



hi 



' /4 ' 



+ 0(e), 
(7.15) 

(7.16) 



In order to be able to express (IT. 12f) in terms of the one-loop operators (IT.4[) — ( T7T6T) it 
is crucial that terms up to order e 2 are included in the latter. 

The coefficients and &qe D describe the running of the electromagnetic coupling 
above and below the electroweak scale, respectively. The former receives contributions 
from all charged fermions and bosons, whereas the latter receives contributions only from 
light fermions, i.e. all charged leptons and quarks apart from the top quark. 

The couplings that enter the one- and two-loop correction factor^ are renormalized 
at the scale /xr. Instead, as discussed in Sect. 16. 1[ the coupling constants in the Born 
matrix element Ai in (I7.3P and (I7.12p are renormalized at the scale Q, i.e. 



M = M 



9i=ai(Q 2 ) 

e=e(Q 2 ) 



(7.17) 



with 



gUQ 2 ) N = 



e 2 (Q 2 ) N - L - 2 ^ 2 



e 2 (/4) 



06,(1)1 

4vr * e 



47T 



XX - 1 



(7.18) 



Thus, by definition, the Born amplitude A4 is independent of the renormalization scale 
/xr, and the dependence of the one- and two-loop amplitudes on /i R is described by the 
terms (I7.13p . 



15 These are the coupling a in the perturbative expansion (|2.5p - (|2.7p and the couplings g\,gi and e 
that appear in (|7.12p - (|7.14p and enter also l|7.4p - (|7.6p through the dependence of the generators (|2. 13|) 
on the couplings and the mixing parameters c w and s w . 
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The contributions (17.131) originate from combinations of UV and mass singularities. 
We observe that the term proportional to vanishes for = M\y. Instead, the terms 
proportional to bi , b% , and &qed cannot be eliminated through an appropriate choice 
of the renormalization scale. This indicates that such two-loop terms do not originate 
exclusively from the running of the couplings in the one-loop amplitude. 

Combining the Born amplitude with the one- and two-loop NLL corrections we can 
write 



M 



NLL 



M ^sew pZ ^em 



(7.19) 



where we observe a factorization of the symmetric-electroweak contributions, 



NLL 



1 + -^Ff 



+ 



Air) 



2 n 



(7.20) 



the terms resulting from the difference between and Mi, 



F z N i^ L 1 + — AF Z , 
4ir 



(7.21) 



and the electromagnetic terms resulting from the mass gap between the photon and the 
W boson, 



NLL 



—AFf 



Aix) 



- (AFf 



(7.22) 



We also observe that the symmetric-electroweak and electromagnetic terms are consistent 
with the exponentiated expressions 



NLL 



cxp 



a 



psew I e \ f~* 



NLL 



exp 



47T 

a 



sew 
2 



AF cm + ( — V 
4tt 1 V4vr/ 



AG^ m 



(7.23) 



In particular, these two contributions exponentiate separately. This double-exponentia- 
ting structure is indicated by the ordering of the one-loop operators F^ ew and AFf 11 in 
the interference term Ff ew AFf m in our result ( 17. 12ft . It is important to realize that the 
commutator of these two operators yields a non-vanishing NLL two-loop contribution. 
This means that the double-exponentiated structure of the result is not equivalent to a 
simple exponentiation, i.e. 



F sew F em ^ exp 



Air 



(F x sew + AFf 11 ) + (^-) 2 (G S 2 CW + AG" 



A-nJ 



(7.24) 



Instead we observe that in NLL approximation the commutator of AF X Z with all operators 
in (I7.19p - (l7.22p vanishes, and also [AFf] 2 does not contribute. Thus, we have 



psew pZ pern NLL pZ psew pern NLL psew pern pZ 



(7.25) 
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and, in principle, F z can be written in exponentiated form, 



P Z NLL / C^t a t-iZ 

r = exp — Ai 1 , 

V47T 



(7.26) 



or absorbed in one of the other exponentials, 



Fsew rnZ NLL 
i< = exp 



pZpem NLL ^ 



47T 



a 



(Ff w + Aif) + 



a- \ 2 
4tt/ 2 



cm 
2 



(7.27) 



The one- and two-loop corrections (I7.3I) - (I7.6I) and (I7.12I) - (I7.13I) contain various com- 
binations of matrices If , which are in general non-commuting and non-diagonal. These 
matrices have to be applied to the Born amplitude Aio according to the definition (12.281) . 
In order to express the results in a form which is more easily applicable to a specific 



process, it is useful to split the integrals I(e, My; —r%j) and A/(e, M\ 



in dZ3D-dZSD 



into an angular-independent part /(e, M v ; Q 2 ) and AJ(e, M v ; Q 2 ), which only involves e 
and L, and an angular-dependent part, which additionally depends on logarithms of r^. 
This permits to eliminate the sum over j for the angular-independent parts of (I7.4I) - (I7.6I) 
using the charge-conservation relation (12.301) . 



(7.28) 



i=l 3=1 



4 = 1 



Moreover, one can easily see that the angular-independent part of (17.41) leads to the 
Casimir operator (12.161) . After these simplifications, all operators that are associated 
with the angular-independent parts can be replaced by the corresponding eigenvalues, 
and the one- and two-loop results can be written as 



m^=m /r+A/r + A/f 



(7.29) 



and 



,~> NLL , , 

Mo = Mr 



[/rf + /r w A/r + - \Mr? + /r aa z + a/? a/? 



em 



(7.30) 



with 



/■sc 

Jl 



A/r 



NLL 



_I ( L 2 + -L 3 e + -L 4 e 2 - 3L - -L 2 e 
2 V 3 4 2 

+ (L + l -L 2 , + ±LV) ^ d + (9(e 3 ), 



Lh 2 



E 



£l fViY ,92 



NLL 
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2e~ 



3e~ 



L 2 



-L a e 



1 



L 4 e 2 



3I + ^ + ^V)E^ 
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- (V 1 + L + \L\ + £ £ l ijm + 0(e 3 



Z NLL 



L + L 2 e + ^ 3 e 2 )/ z £(^c w t 3 - 



sew NLL 
92 — 



Ag™ N = L 



-l 3 - -l r l 2 ) y 

6 2 m I ^ 



9i u {i) (Vi 



_-2 



i=l L w 



5i ^ 



32,(1), 



2 L Q lm 



+ 0(e 3 ) 
+ 0(e), 



+ (^- 3 + + ^V 1 ) 6« D } £ + 0(e) 



(7.31) 



where Z MR = In (/z^/M^), and c i; t 3 , yi, q iy represent the eigenvalues of the operators Cj, 
T 3 , Yi, and Qi, respectively. The only matrix- valued expression in (I7.3ip is the angular- 
dependent part of the symmetric-electroweak contribution ff ew , 



>Q d = EE^ E % ij. 

i=l 3=1 V=A,Z.W ± 



(7.32) 



The two- loop corrections f)7.30p involve terms proportional to /Q d and [/Q d ] 2 . However, 
the latter are of NNLL order and thus negligible in NLL approximation. The combination 
of the matrix (I7.32p with the Born amplitude, 



M)£i d = EE^- E aC" 

i=l 3=1 V=A,Z.W ± 
n n ( 



ipi...tp' i ...ip'....tp„ y y 



= 1 3=1 
3& 



9l ViVj 92.3.3 

r t ,• t-i 



e 2 4 



o2 l J 



+ E AC" 



(pi...tp i ...ip r ..lp n y y 



(7.33) 



requires the evaluation of matrix elements involving SU(2)-transformed external fermions 
ip'i, ip'j, i.e. isospin partners of the fermions ip i7 <pj. 

8 Discussion 

In this section we compare our results presented in Sect. [7] to existing results from the 
literature and apply them to specific processes. 

8.1 Extension of previous results 

In Ref. [27] the one- and two-loop LL and angular-dependent NLL contributions have 
been calculated for arbitrary non-mass-suppressed Standard Model processes, using a 
photon mass and fermion masses for the regularization of soft and collinear singularities, 
respectively. The purely symmetric-electroweak parts of our results, made up of Fl ew or 
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fl™ and its square, confirm the corresponding terms involving 5 sew in Ref. [27], which can 
be seen e.g. by comparing (I7.29p . (I7.30P and (17.311) from our paper with (4.9), (4.10) and 
(4.24) from Ref. [27]. In addition to the existing results we have added all the remaining 
(non-angular-dependent) NLL contributions, including those proportional to ln(M|/M^) 
and the terms involving /3-function coefficients, as well as higher orders in e in the one-loop 
results. We cannot compare the electromagnetic parts contained in AF° m or Af± m and 
<5sem because of the different regularization schemes for the photonic singularities, i.e. for 
the soft and collinear divergences resulting from massless photons. 

In Ref. [28] the complete one- and two-loop LL and NLL contributions for the elec- 
troweak singlet form factors have been derived. We can easily reproduce these form factor 
contributions as a special case of our results, as for only two external fermions no angular- 
dependent logarithms appear and the summation over external legs, using (17.281) . is trivial. 
The functions /, AJ, J, and AJ appearing in (4.69), (4.70), and (4.73) of Ref. [28] cor- 
respond to the equally named functions from this paper for = —Q 2 = s. In the form 
factor case, due to the absence of the angular-dependent term /Q d in (I7.3ip . the one- loop 
operators F[ ew and AF± m commute, so that the two-loop result can be written as a single 
exponential and we have M 2 N = L M H[Fl ew + AF° m + AFf} 2 + G| cw + AG c 2 m }, corre- 
sponding to the form in Ref. [28]. Note that we have renormalized all couplings of the 
loop corrections at the unique scale /i w = /i e = /i R , where the last two lines of (4.73) in 
Ref. [28] vanish. 

8.2 Comparison to Catani's formula in QCD 

The structure of our results for electroweak logarithmic corrections is similar to the 
singular structure of scattering amplitudes in QCD. In Ref. [35] the singular part of a 
QCD one-loop amplitude, i.e. the pole part with terms e~ 2 and e _1 , is obtained from 
the Born amplitude by applying an operator I^\e,fi 2 ; {p}) which, for only fermions and 
antifermions as external particles, can be expressed through Eqs. (12)-(15) of Ref. [35] 
as 

asW) *»(«, Ai M) = L ^ (4) EST. ' Tj/feOi -r„), (8.1) 



2tt v ' rft ' L " J/ 4tt V 2 



where J(e, 0; — r^) is the function defined in (17.71) . and as, e represents the strong coupling 
renormalized at /i R and rescaled as in ( 12.71) . The product Tj ■ Tj of the colour charge 
operators corresponds to J2v=a,z,w ± If Ij f° r electroweak interactions, so that (18.11) has 
exactly the form of F( ew ( I7.4p for M\v = 0, including the angular-dependent logarithms. 
Note that the factor (Q 2 / /i R ) e from the difference in the definitions of as [35] and a t ( 12.71) . 



«s(/4) = «s,, — , (8.2) 




is multiplied by a factor (— ^/rij) 6 from I^\e, /i R ; {p}), producing the correct angular- 
dependent logarithms Uj = ln(— r^/Q 2 ) contained in J(e, 0; —rij). The e-dependent pre- 
factor in the definition of J^(e, /i R ; {p}), e~ e ^^ /T(l — e) = 1 + 0(e 2 ), is irrelevant in 
NLL accuracy, and the same holds for the corresponding factors in the two- loop result. 
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In two loops, the operator acting on the Born amplitude is given by Catani's formula, 
Eqs. (18)-(21) of Ref. [35], in NLL accuracy as 

~ f/ (1) (e,/4; M)f + ^ [i (1) (2e,/i 2 R ; M) - / (1) ( e ,/4; M)l 1 ■ (8.3) 




The first term in the curly brackets originates from the combination of J*- 1 ** (e, /x R ; {p}) 
applied to the one-loop amplitude and — |[J^(e, /x^; {p})] 2 acting on the Born amplitude 
in Ref. [35]. It can be identified with the term |[-Fi sew ] 2 from our two-loop result (17.121) . 
Using colour conservation, the second term gives 

(^) a ^[l«(2e,A4;W)-/ (1) (e,A{p})] 

^ L (^) 2 ^ -WTU(e,0,^), (8.4) 



WW ,u 2 

where J(e, 0, /x^) is the function defined in ( 17. 14ft . In the electroweak model, the expression 
An(3 T 2 in (E3D corresponds directly to the term [b^ g 2 {Yi/2) 2 + b^ ] gld}/ e 2 in (173311 . 
Therefore (18.41) can be identified with the contribution of G| cw for M\y = 0. 

The symmetric-electroweak part of our results can thus be obtained from Catani's 
formula by an obvious replacement of gauge-group quantities together with a simple 
substitution of massless one-loop integrals by massive ones. The remaining parts of our 
results, which are due to the differences between the masses of the photon, Z boson, and 
W boson and which cannot be inferred from Catani's formula, can be expressed as simple 
combinations of the same one-loop integrals. 

8.3 Comparison to electroweak resummation results 

In Refs. [18-21] a resummation of electroweak one-loop results has been proposed, 
including all LL and NLL corrections apart from the ln(M|/M^) terms. 



The non-angular-dependent (9(e°)-terms of /f™ are found in Ref. [18] in the sum 



n f 
k=l 

of Eq. (48) [Eq. (49) in the hep-ph version] for the contribution of external fermionic lines 
above the weak scale. These expressions are extended in Eq. (48) of Ref. [20] in order to 
include the /3-function terms of gf™. Here care must be taken to use the resummed one- 
loop running of the couplings, i.e. a(s) = a(M 2 )/(l + c In s/M 2 ) with c = a(M 2 )(3 Q /TT, 
in order to correctly arrive at 

1, s (. a(M 2 ) \ 1, a(M 2 ) 1 2 s c. 3 s 2 . 

- In — In \ / - 1 + — In \ / = - In 2 — - - In 3 — + 0(c 2 ), (8.5) 

c M 2 \ a(s) J c 2 a(s) 2 M 2 6 M 2 K J ' y ' 

and equivalently for a — > a', c — > c', reproducing the coefficient of the /^-function terms 
in sf w (ESQ) for /x R = M w . 

The angular-dependent NLL corrections are treated in Ref. [21]. The last line of 
Eq. (13) there [Eq. (11) in the hep-ph version ] matches the JCf-tenn of (17311 and (17321) 
in O(e ) if one takes all gauge-boson masses to be my a = Mw- The symmetric-electroweak 
parts of the one- and two- loop amplitudes in Ref. [21] are thus in agreement with our 
results. 
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To summarize, our explicit one- and two-loop results (17.291) . (I7.30P and (17.311) confirm 
the symmetric-electroweak parts of the resummed ansatz. The electromagnetic contribu- 
tions, i.e. the contributions from below the weak scale, cannot be compared due to the 
different regularization schemes for the photonic singularities. However, the factorization 
of the electromagnetic contributions from the symmetric-electroweak ones and the fact 
that the former can be expressed in terms of one-loop QED corrections are in agreement 
with the approach proposed in Refs. [18-21]. 

8.4 Four-fermion scattering processes 

We now apply our results to massless four-fermion processes 

<Pi(px) -> <Pa(-pa) <Pi(~P4), (8-6) 

where each of the may be a massless fermion, ipi = or antifermion, ifi = f£j, with 
the notations from Sect. [21 provided that the number of fermions and antifermions in the 
initial and final state is equal. The scattering amplitudes for the processes (18.61) follow 
directly from our results for the generic n — > process (12.11) by crossing symmetry. The 
Mandelstam invariants are given by s = r 12 = r 34 , t = r 13 = r 24 , and u = r 14 = r 23 with 

Tij = (Pi+Pj) 2 - 

The following discussion focuses on s-channel processes of the form 

K ft ^> tiff (8-7) 

with one fermion and one antifermion in the initial as well as in the final state, where /* 
and /*/ are neither identical nor isospin partners of each other, and the same holds for f* 
and Therefore the external fermion lines are always connected between /" and fp in 
the initial state and between f£, and /5 in the final state. The number of independent 
chiralities Ki is thus restricted to two, k and k' , and the particle pairs in the initial and 
final state must either be antiparticles of each other or antiparticles of the mutual isospin 
partners. The first case corresponds to neutral-current four-fermion scattering, which is 
treated in Sect. 18.4.11 The second case refers to charged-current four-fermion scattering 
and is treated in Sect. 18.4.21 

The scattering amplitudes of all other processes (18. 6p can be obtained from (18.71) 
by crossing symmetry and additive combinations of amplitudes. For instance, the full 
four-quark amplitude for uu — > dd is given by the sum of the s-channel neutral-current 
amplitude and the t-channel charged-current amplitude, where the latter follows from 
the result of the s-channel charged-current amplitude for ud — > ud by exchanging p 2 
and P3, and thus s and t. The neutral- and charged-current results in Sects. 18.4.11 and 
18.4.21 together are thus sufficient to construct the scattering amplitudes for all massless 
four-fermion processes (18.61) . 

8.4-1 Neutral- current four-fermion scattering 

This section deals with the s-channel neutral-current four-fermion processes 
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where a fermion-antifermion pair annihilates and produces another fermion-antifermion 
pair. The electromagnetic charge quantum numbers of the external particles are given 
by <?/ = qf S = -qj 5 and q r = q fK > = -qjj, the hypercharges by y f = y fs = -yj 5 and 
UP — Urx.' — —UfK'i the isospin components by t 3 f = t% = —t% and t% = t 3 , = —t% K ,, 

and the isospin by tf = \t 3 \, tji = \t 3 ,\. These electroweak quantum numbers depend on 
the flavour and the chirality of the fermions. 
The Born amplitude reads 



s 

in the high-energy limit, where 

-4nc = v(p 2 , k)7 m m(pi, k) u(-p 3 , n')l^v(-p 4 , k') 



19) 



.10) 



represents the spinor structure of the incoming and outgoing massless fermions with chi- 
ralities k and k' and 



Co,nc — g\{Q 



VfVf 



+ g 2 2 {Q 2 )t 3 t 3 ,. 



(8.11) 



As indicated, the couplings gi{Q 2 ) in the Born amplitude are renormalized at the scale Q, 
whereas the additional couplings and mixing angles in the loop corrections below are 
renormalized at the scale /xr. 

We write the neutral-current amplitude in the form 



M 



NLL 



NC 



KA psew rpL pern 

• /vl 0,NC- r NC r NC r NC 



(8.12) 



according to (17.191) . Applying the non-diagonal operator _p sew f|7.20p to the Born ampli- 
tude, we find 



-^0,NC F NC ~^ NC l ^0,NC + ^ 



3 4 



-L 3 e 2 



x C ,nc C% c + ( L + + h 3 ,^ Q d NC + 0(e 3 ) 
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AttJ 
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-L 4 - 3L 3 j C 0i nc (C%nc) - L 3 C± d NC C^c + Co,nc #2?nc + ^( e ) 



(8.13) 



where 



(8.14) 



results from the non-angular-dependent contributions to ff ew in (I7.3ip . whereas the ap- 
plication of the angular-dependent operator /Q d on the Born amplitude yields 



C\ 



ad 
NC 



Co. 



NC 



,u, D +in (i 



.3.3 
t f t f 



(8.15) 
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The last missing part in (18.131) . 



n 

i/2,NC 



NLL 



+ C(e), (8.16) 



results directly from (17.31!) . where the values for the /3-function coefficients &i and 6^ 
in the electroweak Standard Model are given in ( 1C.7I) . 

The symmetric-electroweak result (I8.13P is multiplied with the diagonal factors 



F z nll 1 a e 2 



,NC 



and 



NC 



^ 2 (A/r 



47T/ 



em 
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em 
2,NC 
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The factors F^ c and F^™ follow from (17.3 ip . 
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(8.19) 
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L 1 
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+ (^ 3 + 2Le~ 2 + L 2 e" 1 ) 6« D } (g 2 . + q}) + O(e) 



(8.20) 



(8.21) 



where the values for the /^-function coefficients and &qed m the electroweak Standard 
Model are given in (107] and fl08l) . 

Our result can be compared to Refs. [22, 23, 31], where the logarithmic two-loop contri- 
butions to neutral-current four-fermion scattering have been determined by resummation 
techniques. These papers use the renormalization scale /xr, = Mw for the couplings in 
the loop corrections, so that we have to set l m = ln^^/M^) = in (18.161) and (18.211) . 
The large electroweak logarithms in Refs. [22,23,31] are defined with the choice Q 2 = —s, 
such that ln(-s/Q 2 ) = in <KT5\i and (IQjjl . 

In Refs. [22, 23, 31] the photonic singularities are regularized with a finite photon mass. 
Therefore we cannot directly compare their electromagnetic corrections with our dimen- 
sionally regularized result. But Refs. [22,23,31] define a finite scattering amplitude by 
factorizing the complete QED corrections, 



A// fin 
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NC 



QED 

NC 



(8.22) 
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The factor U^q D represents the full QED corrections from photons coupling exclusively 
to fermions, with the photonic singularities regularized in the same way as in the ampli- 



tude -M NC . The soft and collinear divergences resulting from massless photons cancel in 
the ratio (KT22I) . and M^ n c is independent of the regularization scheme for the photonic 



singularities. Using our dimensional regularization scheme and /ir 
above, U§q D is given by 
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In contrast to i 7 ^™ 



U^q U , according to its definition in Refs. [22,23,31], contains 



3 J ^ 3 

(JEH rrQED 

the complete photon contribution without the subtractions at a photon mass equal to M\y. 
This is why ( I8.24p and (I8.25|) differ from (18.201) and (I8.2ip by finite logarithmic terms. 
The expressions /i nc an d S^nc can De obtained from AF* m (17. 5p and AGf m (I7.13P by 
replacing A/ — > / and A J — ► J and adding the term proportional to 6^op i n (I8.25p . where 
6to P is defined in ( 1C.9I) and corresponds to the top-quark contribution to the electromag- 
netic /3-function. In our result, this latter term is implicitly contained in g^c ([SUB]) and, 
by construction, cancels in the subtracted expression Ag^c (18-211) . The factor U$q B in 
(I8.23P is valid for m t ~ Mw, whereas the corresponding factor in Refs. [23,31] is only 
valid for a massless top quark. However, the finite amplitude .M^c * s independent of the 
top-quark mass in NLL approximation. 

The one- and two- loop contributions to .M^c (I8.22p can be expressed as follows: 



fin 
NC 



NLL 



-A 



NC 



-Co, 



NC 



L 1 



3L CT a -A£ 



em 
NC 



fS 



QED 
NC 



+ L Ci 5 NC + Co,NC ^/l Z NC 



0(e), 



M 



fin 

2,NC 




r 2 q j \ sisew \ rem i fQED 

Lj — 6±j I 1NC — ^AJ ljNC + J^NC 



_ {L C^NC + Co,NC ^A Z Nc) ^LNC _ A/l.NC + /^nc) 

+ Co,nc (^nc + A^c - ^2 Q nc) } + 0(e). 



(8.26) 



The one-loop result reproduces Eq. (50) from Ref. [22], and the two- loop result agrees 
with Eqs. (51), (52) of Ref. [22] and Eqs. (51), (52), (54) of Ref. [23] for N g = 3 families 
of leptons and quarks. Note that in these papers (a/47r)(g|/e 2 ) is used as the parameter 
of the perturbative expansion. The corrections proportional to A/ 1 Z NC , which account 
for the mass difference of the heavy electroweak gauge bosons, M z ^ M w , have to be 
compared to Eq. (62) of Ref. [31], where the first order of an expansion in the parameter 
5 M = 4 = 1 - M^/Ml is presented. Using l z = ln(M|/M^) = s 2 w + 0(4), Ref. [31] 
gives indeed the first order in of our result, provided that we set s w = and c w = 1 
in A/f NC ( 18.191) . thus neglecting higher orders in 4- 
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8.4-2 Charg ed- current four-fermion scattering 

In order to complete our predictions for massless four-fermion scattering, we apply our 
results to the s-channel charged-current processes 

tifi^tifi, (8.27) 

where the fermions and are the isospin partners of /f 1 and f~b, respectively. The 
hypercharge quantum numbers of the external particles are given by yj = y^ = —yjv 
and VP = y f h = —VfL and the isospin components by t 3 = t% = t% = t% = t%. All 

'a' 1 p 1 ' Jo J p J a i J i 

external fermions have to be left-handed, so t = \t 3 \ = 1/2. 

In the high-energy approximation, the Born amplitude reads 



AVcc = - s Acc ^ (8.28) 



with the spinor structure 

A CC = Ufa, L)7 m m(pi, L) u(-p 3 , L)7 At u(-p 4 , L). (8.29) 
As in the previous section, the amplitude is written in the form 

M cc N ^ L M 0tCC F c e c w F c z c Fee • (8-30) 
Applying the non-diagonal operator _F sew ( I7.20j) to the Born amplitude, we find 



u psewNLLl, f gj{Q 2 ) Me 

•^o.cc^cc - -^cc^ — 2 — + 



1} + 2 L 3 £ + I L 4 £ 2 _ 3L _ 3 _ 1 L 3 e 2 

3 4 2 2 

gl{Q 2 ) 



+ ■ ■ 



a, N 2 



sew \ r 3 /^yad ^isew 
-L - 6L I ^i,ccj _ ^ °1,CC°1,CCH 2 92,CC + U \ e j 



(8.31) 



with 



e 2 U 2 J 1,00 



c . .M|!)[ 41n(? )| T _ 2 ( ln Q )+lne 

+ 21n (2) s f ( Q. ) SffltM, 

«c ^ (ji 3 " '„i 2 ) (f| ^ + | f| 4") • (8.32) 



The diagonal factors 



F^c N = 1 + ^ AA Z CC (8.33) 
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and 



TTiem NLL 1 C^e » „ cm 

are expressed through 



a^y 

4tJ 



^(AA C ,cc) 2 + A^cc 



(8.34) 
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CZcc- (8-36) 



9 Conclusion 



We have studied the one- and two-loop virtual electroweak corrections to arbitrary pro- 
cesses with external massless fermions in the Standard Model. In the high-energy region, 
where all kinematical invariants are at an energy scale Q that is large compared to the 
electroweak gauge-boson masses, we have calculated mass singularities in D = 4 — 2e di- 
mensions taking into account all leading logarithmic (LL) and next-to-leading logarithmic 
(NLL) contributions. This approximation includes all combinations a l e~ k \n j ~ k (Q 2 /M^) 
of mass-singular logarithms and 1/e poles with j = 2/, 2/ — 1 and 2/— 4 < k < j. All masses 
of the heavy particles have been assumed to be of the same order M\y ~ Mi ~ Mh ~ m t 
but not equal, and all light fermions have been assumed to be massless. 

The calculation has been performed in the complete spontaneously-broken electroweak 
Standard Model using the 't Hooft-Feynman gauge. All contributions have been split into 
those that factorize the lowest-order matrix element and non-factorizable parts. The non- 
factorizable parts have been shown to vanish in NLL approximation owing to collinear 
Ward identities. All factorizable contributions have been evaluated using a suitable soft- 
collinear approximation and minimal subtraction of the ultraviolet singularities. Explicit 
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results have been given for all contributing factorizable Feynman diagrams. The two- 
loop integrals have been solved by two independent methods in NLL approximation. 
One makes use of sector decomposition to isolate the mass singularities, the other uses 
the strategy of regions. The fermionic wave functions are renormalized on shell, and 
coupling-constant renormalization is performed in the MS scheme, but can be generalized 
easily. 

In order to isolate the effects resulting from the mass gaps between the photon, the 
W boson, and the Z boson, all contributions have been split into parts corresponding to 
Ma = M z = M w an d remaining subtracted parts associated with the massless photon 
and the Z boson. By combining the results of all diagrams we found that the electroweak 
mass singularities assume a form that is analogous to the singular structure of scattering 
amplitudes in massless QCD. The sum of the two- loop leading and next-to-leading loga- 
rithms is composed of terms that can be written as the second-order terms of exponentials 
of the one-loop contribution plus additional NLL contributions that are proportional to 
the one-loop /5-function coefficients. All terms can be cast into a product of three expo- 
nentials. The first inner exponential contains the part of the corrections corresponding 
to Ma = M z = M w , i.e. the SU(2) x U(l) symmetric part. The second exponential 
contains the part originating from the mass gap between the Z boson and the W boson 
and contains only terms involving ln(M|/M w ). The third outer exponential summarizes 
the contributions that originate from the mass gap between the photon and the W boson 
and corresponds to the QED corrections subtracted by the corresponding corrections with 
Ma = Mw- While the second exponential commutes with the other two and does in fact 
get no second-order contribution in NLL approximation, the first and the last exponential 
do not commute. 

If one neglects the NLL contributions proportional to m(M|/M w ), our result confirms 
the resummation prescriptions that have been proposed in the literature. These prescrip- 
tions are based on the assumption that, in the high-energy limit, the electroweak theory 
can be described by a symmetric, unmixed SU(2) x U(l) theory, where all electroweak 
gauge bosons have mass M w , matched with QED at the electroweak scale. Indeed, apart 
from the terms involving ln(M|/M w ), in the final result we observe a cancellation of 
all effects associated with symmetry breaking, i.e. gauge-boson mixing, the gap between 
Mz and Mw, and couplings proportional to the vacuum expectation value. This simple 
behaviour of the two-loop NLL corrections is ensured by subtle cancellations of mass sin- 
gularities from different diagrams, where the details of spontaneous symmetry breaking 
cannot be neglected but have to be taken into account properly. 

In massless fermionic processes, the symmetry-breaking effects are restricted to a small 
subset of diagrams, since the Higgs sector is coupled to massless fermions only via one- 
loop insertions in the gauge-boson self-energies. Using the techniques developed in this 
paper, which are to a large extent process independent, we plan to extend our study of 
two-loop NLL mass singularities to processes involving heavy external particles that are 
directly coupled to the Higgs sector. 

As an application of our results for general n-fermion processes we have presented ex- 
plicit expressions for the case of neutral-current and charged-current 4-fermion reactions. 
In the former case we found agreement with existing predictions obtained with the help 
of resummations prescriptions. 
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Our results are also applicable to reactions that involve massless fermions and (hard) 
gluons, such as 2-jet production at hadron colliders, since gluons do not couple to elec- 
troweak gauge bosons. 



A Loop integrals 

In this appendix, we list the explicit expressions for the Feynman integrals that con- 
tribute to the one- and two-loop diagrams discussed in Sects. 15.11 and 15.21 In order to 
keep our expressions as compact as possible we define the momenta 

fa = Pi + h, k 2 =Pi + k, k 3 =pi + h + l 2 , 

qi = Pj - / 1; q 2 = Pj - l 2j q 3 = Pj -lx- l 2 , l 3 = -h - h-, 

n = Pk- h, r 2 =p k - h, r 3 = p k -h + h, U = h~ h- (A.l) 
For massive and massless propagators we use the notation 

P{q,m) =g 2 -m 2 + i0, P(q) = q 2 + iO, (A.2) 
and for triple gauge-boson couplings we write 

r^^{h, i 2 , h) = g ^(h - z 2 r + g^ 3 (h - hY 1 + g^ih - hY 2 - (A.3) 

The normalization factors occurring in (\2.7\i are absorbed into the integration measure 

,2 \ D/2-2 , D 
l _D_\ A-D a 

and for the projection introduced in f)3.24p we use the shorthand 



H, = (4.)= ( A - D ^o = A (e-QV) 2 " D/2 i%, (A.4) 



n y (r) = — 1l (Tu^j) = -!- Tr (rrfrfj) , (A.5) 

where the second equality holds if T does not involve 75 or LOn^, as it is the case in the 
following equations. With this notation we have 

L>o (mi = [ dfi— 4 ! I 1 / 9 ! x p / v 

J P{h,mi)P{ki)P{qi) 

n 1 \ [ aJ ai -I6(faqi)(faq3) 

D 1 (m 1 ,m 2 ;r ij ) = / dhdl 2 



P(h, m x )P{l 2 , m 2 )P(fa)P(k 3 )P( qi )P(q 3 ) 
-I6(faq 3 )(k 3 q 2 ) 



D 2 (mi,m 2 ;rij) = dhdl 2 

J P(/i, mi)P(Z a , m 2 )P(k 1 )P{k 3 )P{q2)P{q3) 

J P(h,mi)P(h,m9)P(lz,mz)P(ki)P(kz)P 
D^m^m^rij) = J d/~id/ 2 

D 5 ( mi ,m 2 ,r v ) - J ^dl 2 



P(h, m 1 )P(/ 2 , m 2 )P(l 3 , m 3 )P(fa)P(k 3 )P(q 3 ) 
P(h, m x )P{l 2 , m 2 ) [P(^i)] 2 P{fa)P(qi) ' 



■ P{l 1 ,m 1 )P{l 2 ,m 2 )P{fa)P{k 3 )P{k 2 )P{q 1 ) 
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D 6 (mi,m 2 ,m 3 ,m 4 ;ry) = f dZ~idZ~ 2 — - — f* 1 9l ^t 



P(h, mi)P(Z 2 , m 2 )P(l 3 , m 3 )P(h, m A )P{k x )P{ qi ) 
x [V m (h, l 2 , Z 3 )r™ 3 (Zi, l 2 , h) + 2Z 2w Zf ] , 

-Ak iqi 



D 7 (m 1 ,m 2 ,m 3 ;r ij ) = / dZidZ~ 2 

J P(l 1 ,m 1 )P(l 2 ,m 2 )P(l 1 ,m 3 )P(k 1 )P(q 1 

D 8 (m u m 2 ,m ?n m^r i j) = J dhdl 2 ' //J 
D 9 (m 1 ,m2,'m 3 ,m4;r ij ) = J dhdi 2 



P(h, m 1 )P{l 2 , m 2 )P(Z 3 , m 3 )P{l l ,m i )P{k 1 )P{q 1 ) 



' P(Z 1; m 1 )P(Z 2 , m 2 )P(Z 3 , m 3 )P(h,m A )P(ki)P( qi ) 
D 10 (mi,m 2 ,m 3 ;r ij ) = D 7 (mi,m 2 ,m 3 ;rij), 
D (m m m m-r )- fdldl ^fgf Tr ( 7w / 27/M / 3 ) 

£ ) ii,mK,m2,m3,m 4 ;r ij ) = — 4D 8 (mi, m 2 , m 3 , m 4 ; r y), 

-16(Zcigi)(fc 3 r 2 ) 



D 12 {m u m 2 ]r ik ) = J dhdl 2 
Di 3 (m 1 ,rn 2 ,rri3-,rij,rik,rjk) = J dZ~idf 2 



P(Z 1; m 4 )P(Z 2 , m 2 )P(h)P(k 3 )P( qi )P(r 2 ) ' 
8^ 1 g 2 V 3 3 r MlM2M3 (— Zi, Z 2 , Z 4 ) 



' P(Zx, m 1 )P(Z 2 , m 2 )P(Z 4 , m 3 )P(h)P(q 2 )P(r 3 ) 
D u (mi,m 2 ; r y , r fe /) = _D ( m i; rij)D (m 2 ; r H ). (A. 6) 



B Relations between loop integrals in NLL approximation 

In the following we list relations between one- and two-loop integrals defined in App. lAl 
These relations are valid after subtraction of the UV singularities and in NLL approxima- 
tion. They have been obtained from the explicit results listed in in Sects. 15.11 and 15.21 and 
are employed in App. [E]in order to simplify the sum over all NLL two-loop contributions. 

As in Sect. 15.21 also in the following the symbols are used to denote generic mass 
parameters, which can assume the values m>i = Mw, Mz,m t , Mh or rrii = 0, and the 
symbols Mi are used to denote non-zero masses, i.e. Mi = My., in- . Mr. 

Combinations of the 2-leg ladder integrals in (15. 5p and (15.81) can be expressed as 
products of one-loop integrals using 

D 1 (m 1 ,m 2 ;r ij ) + D 2 (m u m 2 ; r„) + D l {m 2 ,m 1 ]r ij ) + P 2 (m 2 , m X ] r^) 
N = L Dofa^ri^D^m^nj) 

N i^ L P (M W ; r ij )P (M w ; nj ) + AP (m i; rij )D (M w ; r y ) 

+ D (M W ; r ij )ADQ{m 2 ] r {j ) + AD {m X ] r i:j )AD (m 2 ; r {j ), (B.l) 

where the subtracted functions AD^ have been defined in (15. 4p . For the relations in 
(IB.ip it is crucial that the results for D (m k ; r^) including terms up to order e 2 are used. 
Moreover, Di and P 2 fulfil the relations 

D 1 (M 1 ,m 2 ',r ij ) N = D 1 (M 1 ,M w ;r ij ), 

D 2 (m 1 ,m2;r i j) = P 2 (m 2 , mi; r i5 ). (B.2) 
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The loop integrals corresponding to the non-abelian diagrams involving two external 
legs (15.101) can be replaced as 

D 3 (M 1 , m 2 , m 3 ; r„) N = L ]-D 2 (m 3 , M x \ r tj ) - D 4 (m 3 , M x ; r i5 ) 

- 6D 9 (m 3 ,M 1 ,M 1 ,m 3 ;r ij ), 
AD 3 (M W , M w , mi; ry) N = L -ADi 2 (M w , mi; ry) - AD 4 (mi, M w ; ry) 

- 6AD 9 (mi, M w , M w , m x ; r„), 
AD 3 (mi, M w , M w ; r y ) N = L AD 3 {M W , m u M w ; r y ) + AL>i(M w , m i; ry) 



AD 2 (M W , mi; r y ) + AL> 4 (M W , m x ; ry) 
1 

2" 



- -ADi 2 (M w ,mi;r 4J ), 



AD 3 (Mw,wii,M w ;ry) N = L -^AD 1 {M w ,m l] r ij ) - AD^M w ,m l]rij ). (B.3) 

The first of these relations has been verified and is needed only if at most one of the 
masses m 2 and m 3 is zero. 

The loop functions for the diagrams (15.131) and (I5.16P are equal up to a minus sign 

D 5 (mi, m 2 ; r y ) N = L -D±{m u m 2 ; r y ). (B.4) 

Next we give relations for the loop integrals appearing in the diagrams with self-energy 
insertions in the gauge-boson line, (J5HHD, (153311 . d52SD, dSSHD, and (I5T3D . The 

loop integrals _D 6 , _D 7 , _D 8 , an d -Dn can be related to D 9 as 

D 6 (Afi, M 2 , M 3 , M 4 ; r y ) N = L 10L> 9 (Mi, M 2 , M 3 , M 4 ; r y ), 
Di lj0 (Mi, M 2 , M 3 , M 4 ; rij ) N = L 4D 9 (M ls M 2 , M 3 , M 4 ; ry-), 

AD 6 (0, M w , M w , 0; r tj ) N = 10A£> 9 (0, M w , M w , 0; ry), 
AD lliO (0, M w , M w , 0; r y ) N = L 4AA,(0, M w , M w , 0; r^), 

AD 6 (0, M 2 , M 3 , M 4 ; r y ) N i± L AL> 6 (M 4 , M 2 , M 3 , 0; r y ) 

N = -12 U ^ U ; U: ' AD 9 (0, M w , M w , 0; r tJ ), 
AD 7 (0, M 2 , M 3 ; r y ) N = L A J D 7 (M 3 , M 2 , 0; r, 



M 2 

NLL -3^AD 9 (0,M w ,Af w ,0;r o , : 



'3 

AD 8 (0, M 2 , M 3 , M 4 ; r y ) N = L AD 8 (M 4 , M 2 , M 3 , 0; ry) 

XLL 3 rAL> 9 (0,M w ,M w ,0;r 



AD 9 (0, M 2 , M 3 , M 4 ; r y ) N ± L AD 9 (M 4 , M 2 , M 3 , 0; r 



\ 1 . 1 /. 1 /. (I- ,. 



1/2 _i_ 1/2 

N = 3 2 ^ 2 3 AD 9 (0, M w , M w , 0; r v ), 
AD llfi {0, M 2 , M 3 , M 4 ; r y ) N = L AD lli0 (M 4 , M 2 , M 3 , 0; r -) 
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NLL g Ml + Mt Q . ^ (B 5) 

The loop integral D 9 can be expressed as 

3D 9 (M w ,M w ,M w ,M w ;r -) N = L -J(e,M w ,Q 2 ), 

3AD 9 (0, M w , M w , 0; r -) N = L - [A J(e, 0, Q 2 ) - AJ(e, 0, M- 2 V )J , (B.6) 

and the loop integral Dn as 

3[A£>u,o(0, 0, 0, 0; nj ) - AL> U)0 (0, M w , M w , 0; r y )] N = L -4A J(e, 0, M 2 ,) (B.7) 

in terms of the functions J and A J defined in (17.141) . 

Combinations of 3-leg ladder integrals in (15.381) can be expressed as products of one- 
loop integrals 

D 12 (m l7 m 2 ; nj) + D 12 (m 2 , mi; r ik ) N = L D (m 2 ; r^) D (mi; r ifc ) 
N = L D (M W ; rij)D (M w ; r ik ) + AD (m 2 ; r i:j )D (M w ; r ik ) 

+ D (M W ; rij)AD (mi; r ife ) + AD (m 2 ; rij)AD (mi; r ik ). (B.8) 

For these relations again the terms of 0(e 2 ) in D (m k ; r^) are crucial. Furthermore 

]T sgn(7r( i ,j,A ; ))D 12 (M w ,M w ;r jfc ) N ^ L 0, (B.9) 

w(i,j,k) 

where the sum runs over all permutations 7t(i,j, k) of k, with sign sgn(7r(z, j, k)). This 
identity simply follows from the fact that D 12 is symmetric under the interchange of i and 
k and does not depend on and rj k . 

The loop functions appearing in the non-abelian diagrams involving three external 
legs (15.411) can be replaced as 

D 13 (M 1 ,M 2 ,M 3 ;r ij ) N = L 0, 

2AZ>i 3 (Mi, M 2 , m 3 ; r ij} r ik , r jk ) N = L 2ADi 3 (mi, M 2 , M 3 ; r ik , r jk , r y ) 

N = L 2ADi 3 (Mi, m 2 , M 3 ; r jfc , r y , r ifc ) 

N = L ADi 2 (M w , m 3 ; r,-&) - AL>i 2 (M w , m 3 ; r ik ). 

(B.10) 

C /3-function coefficients 

In this appendix we give relations and explicit expressions for the one-loop /5-function 
coefficients by^y 2 , b 2 , b^\ &qe D , an d &top that have been used in the calculation. For 
more details we refer to Refs. [8,28]. 

The matrix of /5-function coefficients is defined as 

b vk = ^Tr y (/^) - inr,(I*I*) - ~ £ £ Tr*.(/^), (C.l) 

d D 15 * k=R,L 
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where I Vi are the generators defined in (12.131) and Try, Tr$, and Tr^ denote the traces 
in the representations for the gauge bosons, scalars, and fermionic doublets, respectively. 
The sum J^vj, runs over all doublets of leptons and quarks including different colours. The 
traces read more explicitly!^! 

v 3 y i =A,z,w ± 

e 2 Tr*(/^) = e 2 £ I* I* (C.3) 



e 2 Tr* K (J Vl / y2 ) = e 2 £ J^/^. (C.4) 



Multiplying by^ V2 with generators and summing over all gauge bosons yields 



e 2 E = (yV + ^^a, (C.5) 

from which the coefficients corresponding to the weak couplings g\ and #2 can be read off. 
The coefficient corresponding to the electric-charge renormalization is given by 

b^^b^ A = clb? + slb?. (C.6) 

The explicit values in the electroweak Standard Model are (1$ = 1) 

6cf„ 6s£, 3 



w 



The QED /5-function coefficient is determined by the light-fermion contributions only, 

i.e. 



where iV/ represents the colour factor, i.e. N£ = 1 for leptons and N£ = 3 for quarks. 
The top-quark contribution to the electromagnetic /3-function coefficient reads 



e = (c.9) 



D Summing up the one-loop contributions 

In NLL approximation, the contribution of all bare one-loop diagrams to the matrix 
element for a process with n external massless fermions is given by (17.21) . which results 
from the factorizable diagrams (13.51) . Using the explicit results presented in Sect. 15.11 we 
can write 

M\ N = L M [if' scw + Aif' cm + Aif' z l (D.l) 
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Note that our normalization for the trace of the scalar fields differs from the one used in Refs. [8, 28]. 
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with 



-inn 

1 



A^r m = -^EE^/A J D (0;r ij; . 

Z j=l i=i 
n n 



Aif ' Z = -^EE^/A J D (M z ;r u ). (D.2) 

Using the charge-conservation identity (12.301) . the counterterms (16.61) can be cast into the 
form 

1 n n ( - ~1 

^r FN =^oEE^ E /f/^(M w ;Q 2 )+^//AC(0;g 2 )| (D.3) 



i=l j=i Kv=A,Z,W ± 



with 



C(M W ; Q 2 ) ^ L + l Lh + \h\ 2 + 0(e 3 ), 
2 6 

AC(0; Q 2 ) N = L —- — L — h 2 e - \l\ 2 + 0(e 3 ). (D.4) 
e 2 6 

When adding flTXTj) and (|D~3|) we find (|Z3J-flZ2D with 

7(e, M w ; -r -) N = L D (M W ; r tJ ) - C(M W ; Q 2 ), 



A/(e,M z ;-r -) N =Afl (M z ;r 



AJ(e, 0; -r y ) N = L AD (0; r. y ) - AC(0; Q 2 ). (D.5) 
E Summing up the two-loop contributions 

In this appendix, the two-loop results listed in Sect. 15.21 are summed and decomposed 
into reducible contributions, which involve products of the one-loop integrals D (15.21) . 
plus remaining irreducible parts. To this end, we split the integrals according to (15.41) 
and use the relations given in App. [B]as well as the commutation relations f!2.29|) and, in 
particular, the fact that If, If , and J2v=a,z,w ± ^k^k commute with each other. As we 
show, all irreducible contributions cancel apart from those that can be expressed in terms 
of the one-loop functions J (17.141) and /3-function coefficients. To start with, we consider 
four separate subsets and combine these in a later stage. 

Terms related to two external lines not involving gauge-boson self-energies 

We begin by considering the contributions that result from the diagrams where the 
soft-collinear gauge bosons couple to two on-shell external lines and that do not involve 
self-energy contributions to the soft-collinear gauge bosons, i.e. from the diagrams 1, 2, 
3, 4, and 5 of Sect. 15.21 



7 



M 3 / j + M 4 / j + M 5 / j + (i ^ j) . (E.l) 
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These can be summarized as 



V 1 ,V 2 =A,Z,W ± I 



+ I? 2 I*lflfD 2 (M Vl ,M V2 ;r 



{ 9l £ e v ^I^I^I^D 3 (M Vl ,M V2 ,M V3 ;n 

e V 3 =A,Z,W ± 
+ iflflflfD^My^My^) 

+ iflflflfD^My^My^) + j) 



NLL 



-^0 



2 



E iriripipD (M w ; rtf) A>(M W ; r y ) 

V=A,Z,W ± 
V=A,Z,W ± 



3 3 
rZ T A tZ tA 



ijj 



+ lfl?lflfAD (M z ; rij )AD o (0; n> 



+ 



E i^£ A ™ (/f^// 2 + I^f/f 2 ) AD 12 (M W , 0; 

Vi,v 2 =a,z,w± e 



+ E ^\e zv ^ + If If I*) AD 12 (M W , M Z ; r«) 

+ (i j 

+ 6 E l ^ l f Trv(I Vl I V2 )D 9 (M w , M w , M w , M w ; r y ) 

Vi,V2=A,Z,W± 

+ 3 E (^/f+^/^Try^^^ADg^Mw^w^;^)}, (E.2) 

where we have made use of the identities (IPO) . flR2l . flRBl . flR4|) . and (jCT . 

Terms involving gauge-boson self-energy contributions 

The contributions where a soft-collinear gauge boson connects two external lines and 
involves self-energy corrections result from diagrams 6-11 in Sect. 15.21 and read 



M% e = E M? ij . 

m=6 



(E.3) 



They can be summarized as 



M 



2,se 



MM 



O p 2 ^ 

Vi,V 2 ,V3,V4=A,Z,W ± 



if if £ ViV2V 3£ v 4 v 2 v 3 D ^ Mvi ^ Mv ^ My^nj) 
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2 2 

— e v 



4 E lfH?e v ™e v ^ v * (D - l)D 7 (M Vl , M V2 , M n ; r y ) 

x D s {M Vl ,M* 2 ,M V3 ,M VA ; 

r ij) 

\ E E 

r ij) 

~ E E {^,^} 

r ij) 

yi,V4=A,z,iy± 



J <j> <j> J <j> <j>. 



Ei J2 I$? m>« i$ D u , (M Vl ,m i2 ,m i3 , M Vi ; r^) 

* l# i21 #i 3 =it,<i k=R,L s J 



NLL 



-7W 



, n £ I^lf 2 b^ V2 J(e,M w ,Q' 

V 1 ,V 2 =A,Z,W ± 



+ Iflfb^ [A J(e, 0, Q 2 ) - A J(e, 0, M 
+ lflfb$l D A J(e, 0, M^) 

+ E 6/f 1 /f 2 Tr y (/ v ' 1 / v ' 2 ) J D9(Mw,Mw,Mw,Mw 

Vi,y 2 =A,z,VK± 

+ Qlflf Try A J D 9 (0, M w , M w , 0; r i3 ) 

+ 3(/^// + If//) Tr v (I A I z )AD 9 (0, M w , M w , 0; r^j, 



(E.4) 



where we used the relations (jR5j) . flR6|) . flR7l) . (RIB . fl02|) . fl03l) . (jCTl . and (jCTl . For 
the simplification of the diagram (I5.26P in addition (I2.19P was employed. We note that in 
this calculation terms of the form Y,v 2 ,v 3 =a,z,w± iflf e A?293 e ZV2V3 M V2 cancel, where M V2 
results either from the loop integrals or from (I2.19p . 

Terms related to three external lines 

The terms where the soft-collinear gauge bosons couple to three of the n on-shell 
external lines result from diagrams 12 and 13 in Sect. 15.21 and can be written as 



M i3k 



i E ^2 ,llk I + ^\ 



(E.5) 
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where the sum runs over all six permutations n(i,j,k) of external lines i,j,k. These 
contributions yield 



Mf = M) £ £ lf 2 I?Ipl 2 D 12 {M Vl ,M V2 -r ih ) 

W(i,j,k) Vi,V 2 =A,Z,W ± 



— 1- 



■ 92 



£ e v ^lfufu^D 13 (M Vl ,M V2 ,M V3 - r ih r lk , r jk ) 



V U V2,V 3 =A,Z,W± 



M £ o E iM^f^DoiMw-r^DoiMw-n,) 

n(i,j,k) I V 1 ,V 2 =A,Z,W ± 

+ £ 7fj^J// fc A L>o(Mw;ri i )A J D (0;r ifc ) 
+ £ /f/f//4 z Do(M w ;r i ,)A J Do(M z ;r ifc ) 

+ ~lfl?lfl£AD o (0- r -) AD (0; r ifc ) + I?lf-I?l£AD (M z ; r {j ) AD (0; r lk ) 



e 2 
.Pa 1 



+ E ^~e z * v % z lJ*I?AD 12 (M w ,Mz-r i:i ) 



p 2 

where we used ( jRgJ) , (jR9]) . and IjBlOjl . 
Terms from four external lines 



(E.6) 



Finally, we have the contributions where the soft-collinear gauge bosons couple to four 
of the n on-shell external lines, i.e. diagram 14 in Sect. 15.21 



These reduce according to ( 15. 44ft directly to products of one-loop integrals 
Mf l = M £ I^I^I^I^D 14 (M VlJ M V2 ;r iV r kl ) 

V 1 ,V 2 =A,Z,W ± 



(E.7) 



NLL 



M £ lfHjHtH^D {M Vl -T l3 )D Q {M V2 -r kl ). (E.8) 

V 1 ,V 2 =A,Z,W ± 



Complete two-loop correction 

The contributions from the above subsets of diagrams can be combined for an arbitrary 
process involving n on-shell external massless fermions according to (13.81) or (17. lip as 



•A^ = ££ 



8=1 3=1 



•Mgno-se + Ml c ) + £ 



n i 

E oMf 1 



k=l 



1 = 1 



(E.9) 
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As a first step, we show that all terms that cannot be expressed by the one-loop functions 
D and J cancel. The terms involving explicit factors Try(I Vl I V2 ) or Try(I Vl I A ), i.e. the 
last two lines in (1E.2[) and the last three lines in (IE. 41) cancel directly if we use the fact 
that the latter trace is only no n- vanishing for V\ = A, Z. The irreducible terms involving 
explicit e Vl Vi Vb tensors appearing in (IE. 21) and (1E.6I) yield 



n n f I __ __ 

-HEE{ E if >(W? + I?I?I?)AD 13 (My r ,0 ] r tj ) 
i=i 3=1 Vi.Va^z.w* e 

3'A 

+ E i^e^^lflp + lflflf)AD l2 {M^ M Z ; r y ) 

n r n 1 - - 

+ E E i^eAWtflfiPADuiMv, 0; r«) 

fe=1 L Vi,v 2 = J 4,z,w± e 



+ E i^^ Wiy2 /f// 2 4 yi AA 2 (M W) M z ;r, 



(E.10) 



These terms vanish upon using global gauge invariance (I2.30p . 

The complete two- loop correction is thus given by the contributions to (1E.2I) . (IE. 6ft . and 
(IE .80 involving products of -D -functions and the terms involving /3-function coefficients 
and J-functions in ( IE. 4ft . These can be summarized straightforwardly as 



n n n n 



M 



I N = >oEEEE^ E iFlM'^DoiM^r^DoiM^ru) 

4 i=l i=i fe=l i=i I Z V r i,V r 2=A,Z,W ± 

37^ l^k 

+ E /f/] / 4 A ^o(M w ;r lJ )A J Do(0;r H ) 

V=A,Z,W± 



+ -if Ifl£lfAD (0; ry)AL> (0; r H ) + If Ifl£lfAD (M z ; r y )AD (0; r« 



- o-^oEE E inr^j^M^a 2 

i=\ 3=1 ly x y 2 =A,Z,W± 

37^ 



+ Itlfb { l\ AJ(e,0,Q 2 )-AJ(e,0,il^) + If If A J(e, 0, M^) L (E.ll) 



In NLL accuracy, this result can be expressed in terms of the lowest-order matrix 
element .M and the one-loop correction factors (ID. 21) as 



A ~,F NLL . / 

= M ( 



1 r 



F,sew 



| pF,sew^ piF,em | ^>F,sew^^iF,Z 



+ 2 



F,em 



+ Aif^Aif' 6 " 1 + G^' sew + AG 



F,sew 



xF.em 



(E.12) 
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Note that in NLL approximation F^' sew anc [ AFf' 6111 do not commute, while AF^' 7, com- 
mutes with the other terms. Using (12.301) . (lC.5p . and (10.61) . the terms resulting from the 
last two lines of (IE. IIP can be written as 



2/^yF,sew 

C WO 



AG 



F,em 



1 n 

-T 

1 



J(e,M w ,Q 2 



^J2Qi\b { e } [AJ(e,0,g 2 )-AJ(e,0,M^)] +6« D AJ(e,0,M^ 



(E.13) 



Adding the term (I6.10P resulting from parameter renormalization with D replaced 
by / and using the definition of J, (17TH]) . the functions G b 2 ' sew and AG b 2 ' em in (1ET2D get 
replaced by the functions G?| ew and AGJf 11 as defined in (17. 13j) . 

Finally, when combining the wave-function counterterm (16 .7p with the parts of (IE. 121) 
involving products of the functions these terms can be written in the form given in 
( I7.12p . In order to arrive at this result we write Ai\ in the form (16.81) and SZjp in the 
form (16.51) . After arranging the Casimir operators in an appropriate order, we can use 
global gauge invariance (12.301) to transform the wave-function counterterm contributions 
to the form needed. 
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